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Abstract. We continue our study on variable order Arnold-Falk-Winther el- 
ements for 2D elasticity in context of both affine and parametric curvilinear 
elements. We present an /i-stability result for affine elements, and an as- 
ymptotic stability result for curvilinear elements. Both theoretical results are 
confirmed with numerical experiments. 



1. Introduction 

1.1. Elasticity problem. Linear elasticity is a classical subject, and it has been 
studied for a long time. The presented research is motivated with a class of time- 
harmonic problems formulated as follows. Let f2 C R™, n = 2, 3, denote a bounded 
domain occupied by the elastic body. Assume that the boundary of 0, T = dfl has 
been split into two disjoint subsets Li and T2, 

r = TTur;, rxnL 2 = 

Subsets Li,L 2 are assumed to be (relatively) open in L. We seek: 

• displacement vector Ui(x), x £ 0, 

• linearized strain tensor ejj(x), x G f2, and 

• stress tensor (x) , x e 57 

that satisfy the following system of equations and boundary conditions. 

• Cauchy's geometrical relation between the displacement and strain, 

(1.1) ey = + x eil 

• Equations of motion resulting from the principle of linear momentum, 

(1.2) - (Tyj - p{x)u?u i = fi(x), ie!! 

• Symmetry of the stress tensor being a consequence of the principle of an- 
gular momentum, 

1 j — ^3 ^ ' *^ 

• Constitutive equations for linear elasticity, 

aij = E ijk i(x)e k u x e 0, 

• Kinematic boundary conditions, 

Ui = u", x G Li 
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• Traction boundary conditions, 

(TijUj = gi(x), x e T 2 

Here: 

• p is the density of the body, 

• /, are volume forces prescribed within the body, 

• gi are tractions prescribed on T2 part of the boundary, 

• u® are displacements prescribed on Ti part of the boundary, 

• rij is the unit outward normal vector for boundary T, and 

• Eijki is the tensor of elasticities. 

As usual, commas denote partial derivatives, and we use the Einstein's summation 
convention. Vector tj := (XijUj is known as the traction or stress vector. 

The geometrical relations imply the symmetry of the strain tensor. The symme- 
try of strain and stress tensors imply then the usual (minor) symmetry conditions 
for the elasticities, 

(1-3) Eijki — Ejiki — Eijik 

The laws of thermodynamics imply an additional (major) symmetry condition, 

Eijki = Ekuj 

and positive defmiteness condition, 

Eijki^ij^ki > 0, V£y = £jj 

The last condition implies that the constitutive equation can be inverted to repre- 
sent strains in terms of stresses, 

(1-4) €kl — CklijCTij 

where Cf-uj = is known as the compliance tensor and satisfies analogous 

symmetry and positive-definiteness properties. 

For isotropic materials, the elasticities tensor is expressed in terms of two inde- 
pendent Lame's constants fj, and A, 

E^ki = n(5ikSji + SuSjk) + XSijSki 

and the constitutive equation reduces to, 

&ij = 2fieij + Xekk^ij 

where 5y denotes the Kronecker's delta. Inverting, we obtain 

a 1 A 

ea - Aa tJ , ey - —a tJ - -^——^a^ 

Of particular interest is the case of nearly incompressible material corresponding 
to A — > 00. Notice that the norm of the elasticities blows up then to infinity but 
the norm of the compliance tensor remains uniformly bounded. This suggests that 
formulations based on the compliance relation have a chance to remain uniformly 
stable for nearly incompressible materials. Another reason of using formulations 
based on compliance relation is for (visco)elastic vibration problems for structures 
with large material contrast. For a more detailed discussion on these problems, we 
refer to the introduction in [19] . 



MIXED FEM FOR 2D LINEAR ELASTICITY 



3 



1.2. Variational formulations for elasticity. 

Dual-Mixed Formulation. We eliminate the strain tensor and combine geometrical 
relations (jl.ip and constitutive equation (|1.4|) in a single equation, 

CkiijCij = -(ukj +ui,k) 

Classical formulation for elasticity is based on satisfying the equation above in a 
strong sense while relaxing the momentum equations. The idea behind the dual- 
mixed formulation for elasticity is exactly opposite, the equation above is relaxed 
whereas momentum equations (jl.2[) are satisfied in a strong sense. 

The final formulation includes building in kinematic boundary conditions and it 
reads as follows. 
(1.5) 

<7 e H(div, : o-ijUj = g t on T 2 , u G L 2 (f2,V) 

/ Cijkicrkinj + / UiUj,j = / UiTijiij Vr G iJ(div, fi, S) : r^nj = on T 2 
Jn Jn Jr t 

a^^Vi-u? I puiVi = / f l v i Vt)GL 2 (fi,V) 
n Jn Jn 

Above, V — R 3 and § denotes the space of symmetric tensor. L 2 (fi, V) denotes the 
space of square integrable functions with values in V, and if (div, fi, S) denotes the 
space of square-integrable functions with values in §, whose row-wise divergence 
lives in L 2 (f2,V). The traction conditions are satisfied in the sense of traces for 
functions from i?(div, Q) (they live in H~ l / 2 (T)). 

The corresponding static case can be derived formally by considering the so- 
called Hellinger-Reissner functional, and the variational formulation is frequently 
identified as the Hellinger-Reissner variational principle. 

Dual-Mixed Formulation with Weakly Imposed Symmetry. The symmetry condition 
is difficult to enforce on the discrete level. This has led to the idea of relaxing the 
symmetric function and satisfying it in a weaker, integral form. This is obtained by 
introducing tensor- valued test functions q with values in the space of antisymmetric 
tensors IK :— {qij : qij — — qji}, and replacing the symmetry condition with an 
integral condition, 

CTygy=0, VgGi 2 (0,K) 

n 

On the continuous level, the integral condition implies the pointwise condition 
(understood in the L? sense), but on the discrete level, with an appropriate choice 
of spaces, the integral condition does not necessary imply the symmetry condition 
pointwise. 

The extra condition leads to an extra unknown. The derivation of the weak form 
of the constitutive equation has to be revisited. We start by introducing the tensor 
of infinitesimal rotations, 

Pij — 2^. Ui -i ~ U Jy i ) 

Upon eliminating the strain tensor, the constitutive equation in the compliance 
form is now rewritten as, 
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Multiplication with a test function r (now, not necessarily symmetric), integration 
over f2, and integration by parts, leads to a new relaxed version of the equation, 



CijklCTklTij — — 


/ u i T ij,j ~+ 


[ 


l 




Jn 


Jr 


Jn 



We obtain a new variational formulation in the form: 
(1.6) 

f a G jJ(div,ft,M) : (TijUj = g 4 on T 2 , u G L 2 (S},W), p G L 2 (fi,K) 

Ci-juOkiTij + Uinjj +PijT i: j = I u^TijUj Vt g i?(div, Q,M) : ryn^ = on T 2 
(JijjVi-uj 2 J puiVi = / Vve£ 2 (f2,V) 

t, ,.'/,,. " VgeL 2 {n,K) 



In the above, i?(div,f2, M) denotes the space of tensor- valued square integrable 
fields with a square integrable divergence. The formulation for the static case can 
again be derived formally by looking for a stationary point of the so-called General- 
ized Hellinger-Reissner functional, and it is frequently identified as the Generalized 
Hellinger-Reissner Variational Principle. 

In the analysis presented in this paper, we consider the static case only and, for 
the sake of simplicity we assume that the body is fixed on the entire boundary <9f2, 

i.e. r 2 = 0. 

To connect with the notations used in [B], [2], we denote by A the compliance 
operator. 

A : M 3 eij — ► a ij = CijuSki G M. 

The operator A maps tensors to tensors, is bounded, symmetric and, for piece- 
wise constant material properties, uniformly positive definite. Symmetry proper- 
ties (|1.3I) imply that A maps S into itself. In the isotropic case, A has the form 

Aa = — (a Tr (a) I 

2/i \ 2/x + n\ 

where A (x) , /x (x) are the Lame coefficients. 

We can rewrite then formulation (11.5[) with a more compact, index free notation. 
Find o~ G H (div, Q; S), and u e L 2 (fi; V), satisfying 

(1.7) / (Act : r + divr • u) dx = 0, r G H (div, O; S) , 

Jn 

/ diver -vdx= I f -v dx, v £ L 2 (Q; V) . 
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Similarly, the formulation (|1.6|) becomes to seek a £ H (div, ft; M), u £ L 2 (ft; V), 
and p £ L 2 (ft; K) satisfying 

(1.8) / {Aa -.T + dWr -u + t :p) dx = 0, r £ H (div, ft; M) , 
Jo 

/ diver -vdx = j f -v dx, v £ L 2 (ft; V) , 
Jo Jo 

/ cr : g dx = 0, q £ L 2 (ft; K) . 

Dual-Mixed Formulation with Weakly Imposed Symmetry for Plane Elasticity. In 
two space dimensions, the skew-symmetric tensors involve a single non-zero com- 
ponent p, 

Us) 

The formulation (|1.8j) reduces to seek er G H (div, ft; M), u £ L 2 (ft;V), and p £ 
L 2 (ft) satisfying 

(1.9) / (Act : r + divr • u - S x t p) dx = 0, re H (div, ft; M), 
Jo 

/ diva-vdx= / f-vdx, v £ L 2 (il;Y), 
Jo Jo 



/ Siaqdx = 0, <? 6 L 2 (ft). 
Jo 



p2x2 



where operator 5*i maps a real 2x2 matrix to a real number. For any a £ 

(1.10) S\(J = CT12 - CT 2 1. 

1.3. Review of the existing work. Besides the Hcllinger-Reissner and general- 
ized Hellinger-Reissner principles discussed in the previous sections, there are other 
variational formulations for elasticity, see [18]. 

Stable finite element discretizations based on (|1.7j) are difficult to construct. A 
detailed description of related challenges can be founded in [M] . A stable discretiza- 
tion based on the strong enforcement of symmetry condition was developed in [2] . 
For the lowest order element, the number of degrees of freedom for stress tensor in 
3D is 162. 

Meanwhile, a number of authors have developed approximation schemes based on 
p| . See [51 El HH [U El HI [HI HH HU EH! EH Eg EH] . For a brief description of these 
methods, we refer to the introduction in |6]. Recently, Cockburn, Gopalakrishnan 
and Guzman have developed a new mixed method for linear elasticity using a 
hybridized formulation based on (II .8[) in [10] . 

The work presented in this paper is based on mixed finite element methods 
developed by Arnold, Falk and Winther for the formulation (|1.8|) in [3 HU |6]. We 
restrict our analysis to plane linear elasticity problem only. 

At the first glance, a generalization to elements of variable order seems to be 
straightforward. And this is indeed the case from the implementation point of 
view. The stability analysis, however, is much more difficult. The canonical pro- 
jection operators defined in [5] do not commute with exterior derivative operators 
for variable order elements, a property essential in the proof of discrete stability. 
In [19], we resolved the problem by invoking projection based (PB) interpolation 
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operators. Unfortunately, the PB operators do not commute with the operator 
S' n _2) invalidating construction of operator Wh introduced in the stability analysis 
presented in [5] We resolved the problem in [19] by designing a new operator Wh- 
However, in our first contribution we were only able to prove the well-definedness 
of Wh for n — 2 with orders varying from to 3. In this paper, we propose PB 
operators of new type and a new operator Wh that enable the ^.-stability analysis 
for 2D elements of variable or Additionally, we analyze curvilinear elements 
and establish conditions for their asymptotic ^-stability. 

1.4. Scope of the paper. An outline of the paper is as follows. Section 2 in- 
troduces notations. Section 3 reviews definitions of finite element spaces on affinc 
meshes. In Section 4, we return to the mixed formulation for plane elasticity with 
weakly imposed symmetry, and recall Brezzi's conditions for stability. In Section 
5, we establish necessary results for proving the stability for affine meshes. We 
design new PB operators and a new operator Wh- In Section 6 we prove that the 
Brezzi conditions are satisfied for affine meshes. The main result for affine meshes 
is Theorem 16. 31 In Section 7, we discuss curvilinear meshes and geometry assump- 
tions. Section 8 reviews definitions of finite element spaces on a reference triangle, 
on a (physical) curved triangle, and on curved meshes. In Section 9, we establish 
necessary results for proving the asymptotic stability for curvilinear meshes. This 
part of work is much more technical than that of Section 5 because of the curvi- 
linear meshes. We explain briefly the substantial difference between the analysis 
for curvilinear meshes and that for affine meshes at the beginning of Section |9] 
In Section 10 we prove that the Brezzi conditions are satisfied asymptotically for 
curvilinear meshes, i.e. for meshes that are sufficiently fine. We conclude the pa- 
per with numerical experiments that illustrate and verify the presented theoretical 
results. 

2. Notations 

We denote by T an arbitrary triangle in R 2 . And let T be the reference triangle. 
We denote the set of subsimplexes of dimension k of T by Afe(T), k = 0, 1, 2, and 
the set of all subsimplexes of T by A(T). Ao(T) consists of all vertices of T, Ai(T) 
consists of all edges of T, and A 2 (T) = {T}. 

For U, a bounded open subset in K™, we define: 

C k (U) = {u£ C k {U) : D a is uniformly continuous on U,V\a\ < k}. 

We also define 

C M (£7) = {u e C k (U) : D a is Lipschitz on f/,V|a| = k}. 

M will denote the space of 2 x 2 real matrices. For any vector space X, we denote 
by L 2 (fi;X) the space of square-integrable vector fields on f2 with values in X. 
In the paper, X will be R, R 2 , or M. When X = R, we will write L 2 (Q.). The 
corresponding norms will be denoted with the same symbol L 2 (i7;M). The corre- 
sponding Sobolev space of order m, denoted H m (Q;X.), is a subspace of L 2 (f2;X) 
consisting of functions with all partial derivatives of order less than or equal to m 



With an upper bound on the polynomial order. 
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ijm(n). The space 



in L 2 (f2;X). The corresponding norm will be denoted by 
H(div, n- M) is defined by 

ff(div,0;M) = {a e L 2 (fl;M) : diver G L 2 (ft;R 2 )}, 

where divergence of a matrix field is the vector field obtained by applying operator 
div row- wise, i.e., 



diver = ( 



dan d(Ji2 da 2 i 



dxi 8x2 ' dx\ 0x2 
We introduce also a special space 

H(Q) = {cue ff(div, Q) : uj\ dn G L 2 (<9fJ; M 2 )} 

with the norm, 

IMIff(Q) = \\u\\ H (div,n) + lklU 2 (9n),Vu; € F(O). 
For any scalar function u and any vector function u) with values in 

duii duii 

curl « = ( — — , — — — ) 1 , curl uj 

0x2 0x1 



we denote 



dx 1 
9w 2 



Finally, by 



dx 2 ' 
du>2 

dx 2 ' dx\ 

we denote the standard 2-norm for vectors and matrices. 



3. Finite element spaces 



We begin by introducing the relevant finite element spaces on a single triangle 
T. Then, we define the finite element spaces on a whole affine triangular mesh Th 
by "gluing" the finite element spaces on triangles. 

3.1. Finite element spaces on a single triangle. For any r G Z + := {n € Z : 

n > 0}, we introduce 

(3.1) V r (T) := {space of polynomials of order r on T}, 

P r A°(T) = V r A 2 (T) := Vr(T),V r A 1 (T) := [P r (T)} 2 , 
■p r A°(T) = V r {T) := {w e V r A°(T) : w\ e = 0,Ve e Ai(T)}, 
P-A^T) := [7Vi(r)] 2 + (x 1 ,x 2 ) T P r _ 1 (T), 
VrA (T;R 2 ) = V r A 2 (T;R 2 ) := [V r (T)} 2 , 

V r A x {T;R 2 ) := 



: (cril,cri2) T ,(cr 2 l,a22) T € 'P r A 1 (r) | . 



Oil CT12 

021 022 

In [5] [14] , spaces in (|3.1[) are defined in the language of exterior calculus. Here 
we just rewrite them using the language of calculus. We refer to [5] and [TJ] for a 
detailed correspondence before the exterior and classical calculus notations. 



We denote by f a mapping from A(T) to Z + such that if e, / G A(T) and e C / 
then f(e) < r(f). We introduce now formally the FE spaces of variable order. 

Definition 3.1. 

VfA°(T) := {u G P f(T) A°(T) : Ve G A X (T), u| e £ Pf (e )(e)}, 
7>,A 2 (T) := P f(T) A 2 (r) = P r - CT) (T), 
P f A x (T) := G Vf [T) A l (T) : Ve G Ai(T), w • n| e G V m {e)}, 



WEIFENG QIU AND LESZEK DEMKOWICZ 



PfA\T) := {lo e V^A^T) : Ve e Ai(T), u ■ n| e € ^ (e) _i(e)} 
P f A°(T;M 2 ) := [7>fA°(T)] 2 , P f A 2 (T;M 2 ) := [7> f A 2 (T)] 2 , 

0U,<T12) T ,(<721,<T22) T eV f A l {T) 



Cll Cl2 
C21 022 

Here n is the outward normal unit vector along 9T. 

Remark 3.2. According to [S], for any e £ Ai(T), we have 

7> f A (T)| e = P He) (e),VrA\T)\ e • n - P m (e),VfA l (T)\ e • n = P f(o) _i(e). 

Actually, spaces in Definition 13.11 have been given in [T§] with the language of 
exterior calculus. 

3.2. Finite element spaces on an affine triangular mesh. Let Th be an affine 
triangular mesh. We extend the f mapping to a global map defined on A(7/j) with 
values in Z + such that if e C /, then f(e) < f(f). 

Definition 3.3. We put fifc := UreTJ ^ • 

CA (Th) '•= {u E H 1 ^^) '■ u is piece-wise smooth with respect to Th}, 

CA 1 (Th) := {oj 6 H(div, £2/,) : ui is piece-wise smooth with respect to Th}, 

CA 2 (Th) '■= {u £ L 2 (Qh) '■ u is piece-wise smooth with respect to Th}- 

We define 

V f A°(T h ) := {u e CA°(T fe ) : u| T g Pj=A°(T),VT e T ft }, 
PfA 1 ^) := £ CA^Tfe) : u\ T € V f A\T),\/T e 7"fc}> 
V~A l (T h ) ■■= W e CA 1 (7^) : w| T e VfA l {T),VT e 7^}, 
P f A 2 (Tft) := e CA 2 (T h ) : u|r G ^A 2 (T),VT e 7/,}, 
^A (T,;R 2 ) := [Pf A°(7h)] 2 , VfA 2 (Th\ K 2 ) := [P f A 2 (7^)] 2 , 



PfA 1 ^: 



CTll 012 
021 022 



(<T11,<T12) T ,(<721,CT 22 ) T G ^(r^K 2 ) 



Remark 3.4. Spaces V?A°(T h ), VfA l {T h ), V?A l (%), VfA 2 (T h ) coincide with those 
analyzed in [19] . 

Lemma 3.5. 

PfA 1 ^) c V^A^Th) c Vf +1 A\T h ), 
diAfPr+xA l {T h ) c V f A 2 {T h ), curtP f+1 A°(Th) c VfA l {T h ). 
Proof. This is straightforward. □ 

4. Mixed formulation for elasticity with weakly imposed symmetry 



We assume that there is 
f (/) < r n 



€ N such that, for any h > and / £ A (7a), 



We recall the mixed formulation (TQ|) : Find (er, u, p) e /f (div, ft; M) x L 2 (£2; 



L 2 (£2) such that 



(4.1) 



(Aj,t) + (divr,w) - (Sxr,p) = 0, re iJ(div, £2; M), 
(diva, «) = (/,«), DSl 2 (fl;K 2 ), 
<Si<7,g)=0, qeL 2 (n). 
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Here (•, •) is the standard L 2 inner product on ft. This problem is well-posed. See 
[5] and [14] for the proof. 

We consider now a finite element discretization of (|4.ip . For this, we choose 
families of finite-dimensional subspaces 

A£(M) C #(div,f2;M),A 2 (M 2 ) c L 2 (ft;R 2 ), A 2 C L 2 (n), 

indexed by h, and seek a discrete solution (ah,Uh,Ph) & A^(M) x A 2 (R 2 ) x A| such 
that 

(4.2) (A7 h ,r) + (divr > tt fc )-(ftr > pfc)=0, r 6 A^(M), 

(diva h ,v) = (f,v), «eA 2 (R 2 ), 
(Sia h ,q) = 0, qeA 2 h . 
The stability of (14.21) will be ensured by the Brezzi stability conditions: 

(4.3) (SI) ||T|| ff / div n . M ) < ci(At,t) whenever r £ A^(R 2 ) satisfies (divr, v) = 

Vw £ A^R 2 ) and (Sit,?) = Vg 6 Aj[, 

(4.4) (S2) for all nonzero (t>, q) £ A 2 (R 2 ) x A 2 , there exists nonzero 

r G A^(R 2 ) with (divT,u) - (SiT,q) > c 2 ||T|| H(div n . M) (|MlL 2 (n ; R 2 ) + \\q\\ L 2 m ), 

where constants C\ and ci are independent of h. 

For meshes of arbitrary but uniform order, conditions (|4.3[) and (14 .4[) have been 
proved in [5] and [14]. In what follows, we will demonstrate that they are also 
satisfied for (2D) meshes with elements of variable (but limited) order. 

Before presenting our proof, we would like to comment shortly on difficulties 
encountered in proving stability for generalizing AFW elements with variable order. 
As we have shown in Section 3, it is rather straightforward to generalize AFW 
elements to the variable order case. The following commuting diagrams are essential 
in the stability proof from 6 ( 5 , 14 , 

(4.5) 

ff 1 ^;" 2 ) ^» £ 2 (^) H 1 (it; R 2 ) ^ L 2 (n) ff^^R 2 ) ^4 ff^^R 2 ) 

ify- + n2 4- K i n 2 i n° 4. nj;- i 

P r "+i Al (Tfc) ^ V r A 2 (T h ) 7V+iA l (Tfc) ^ V r A 2 (T h ) V r+2 A°(T h ; R 2 ) ^-4"* 7'- +1 A 1 (Th) 

When uniform order r is replaced by variable order f, then the left and middle 
diagrams do not commute if n^' and are natural generalizations of canonical 
projection operators introduced in (6j [5j [14]. A counterexample is given in the 
appendix of |19) . To our best knowledge, the only operators which make these two 
diagrams commute for meshes with variable order, are projection based interpola- 
tion operators introduced in [19] . But then we need to construct an operator 
that makes the right diagram commute if n^ : ~ is the projection based interpola- 
tion operator. We replaced the operator n° with a new operator Wh, constructed 
in [19], but we were able to prove its well-definedness only for < f(T) < 3. In 
the following sections we pursue a different strategy, defining new projection based 
interpolation n^' - ,!!^ operators, and a new operator Wh in such a way that we 
can make all three diagrams commute. 
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5. Preliminaries for the proof of stability 

We begin by recalling our assumptions on the domain and meshes: Q is a poly- 
gon and it is meshed with a family (7h)h of affine triangular meshes satisfying 
assumption of regularity. For any mesh 7~h, mapping f : A(7fc) — > Z + defines a 
locally variable order of discretization that satisfies the minimum rule, f(e) < f (/) 
for e C /, e, / S At- The maximum order is limited, i.e. swp h sup Te7 - ?{T) < oo. 

Definition 5.1. For any T G Tt, we define a linear operator II? T as the standard 
orthogonal projection operator from L 2 (T) onto VfA 2 (T). 

5.1. Projection Based Interpolation onto 7 , f+iA 1 (T). 

Definition 5.2. For any T G Th, we define a linear operator III, ^ T : H l (T\ M 2 ) — ► 
Pf+xA^T) by the relations 

(5.1) f div(nl +liT a;-a;)(x)V(x)dx = V^eP f(T) (T)/M 

(5.2) / (nl +1 T w(x) - w(x)) T curkp(x)dx = V<y9 G P f . {T)+2 {T) 
Jt 

(5.3) J[(T^ +liT u-w)-n]r}(a)da = V77 G P f(e)+ i(e) ,Ve e A X (T) 

Here n is a unit outward normal vector along e. 

The operator III , x T is the Projection Based Interpolation operator onto V? +1 A 1 (T) 
defined in [19] . We immediately have the following two lemmas. 

Lemma 5.3. For any T G Th, and any lu G [iJ^T)] 2 , we have 

II 2 T divoj = c?swIF~ +1 t oj. 

Lemma 5.4. There exists C > swc/i </iai 

||nJ iT w|| ia(r) < CIIwHhicd VT eVe H X {T\ R 2 ). 

5.2. Modified Projection Based Interpolation onto 'P f+1 A 1 (T) and modi- 
fied operator W onto Vf+2-h°(T). In order to prove the stability of the mixed 
FE method, we need to make the left and right diagrams in (I4.5[) commute. In [115] , 
the definition of Projection Based (PB) interpolation operator onto Vf +1 A 1 (T) was 
very similar to ■Pf+iA 1 (T). From the proof of Lemma 10 in [TIJJ, we can see that 
only conditions (|5.ip and (|5.3[) are used to prove the commutativity of the middle 
diagram in (|4.5I) . This implies that we may be able to change condition (|5.2I) for 
the PB interpolation operator onto , P f ^ 1 A 1 (T) in such a way that we can design 
a new operator Wh which makes the right diagram in (I4.5[) commute, and can be 
proved to be well-defined, for an arbitrary order of discretization. 

Definition 5.5. Let f : A(T) — > Z + be a mapping that prescribes the local order 
of discretization and satisfies the minimum rule, i.e. if e, / G A(T) and e C / 
then f(e) < r(/). We put fcf = dimcurlx'P-^^ +1 (T). Let {ff^x), • • • , ff,k,.(x)} be 

a basis of curlx'P-(j.- )+1 (T). Let {gf.i(x), • ■ ■ jgr.fcfCx)} be a linearly independent 
subset of V f( f, 1 (f;R 2 ) such that 

span{g f , 1 (x), • • • ,g fikf (A)} egradxP-^JT) = [7Vt)-# )] 2 
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For t G [0, 1], we define hf,i(x,t) = (1 - i)f>,i(x) + igf,i(x), 1 < i < k? . 

Remark 5.6. It is easy to check that k? = dim'Pfff\_ 1 (T- ) R 2 ) — dimgrad^'P-^j^T). 

Remark 5.7. For any T G 77i, there exists an affine mapping from the reference 
triangle T onto T, which can be written as x = Bt% + b. Here Bt is a constant 
nonsingular matrix in R 2x2 . We denote by f a mapping from A(T) to Z + such that 
if e, / € A(T) and e C / then r(e) < f (/). With the same symbol f we denote the 
corresponding mapping from A(T) to Z + , f(/) := f(/) for any / G A(T), where / 
is the image of / under the affine mapping mentioned above. 

Definition 5.8. (One-parameter family of PB interpolation operators onto Vf, 1 A 1 (T)) 
For any t G [0, 1], we define a linear operator 11*'" f { : ff^TjR 2 ) — > V^ +l K 1 {f) 
by the relations 

(5.4) / div *( n J+ 1 ,f,t c2 ' ~ = G V Hf) (f)/R 

(5.5) JSK+1 f t^(*) ~ ^(x)) T h l (x,t)dx = 1 < t < fc f 

(5.6) f [(nl'- ift 6j-6j)-h]fids = Vf?GP f( e)(g),VeG A^T) 

Definition 5.9. (One-parameter family of PB interpolation operators onto Pfq_j A 1 (T)) 
For t G [0,1], and for any T £ Ti, we define a linear operator II f \7 X T t : ff 1 (T; R 2 ) — >■ 
Vf+^iT) by the relations 

(5.7) f div(ni; 1Tt w-w)(x)^(x(x))rfx = V|eP fm (f)/M 

(5.8) / (n^ 1Tt w(x) -w(x)) T S T T h i (i(x),t)dx = 1 < i < k? 

JT ' ' 

(5.9) f[Q4+i, Tl t u ~ w ) • n]v(s)ds = VjjG? f(e) (e),VeG A^T) 

e 

In the above, x = x(x) denotes the inverse of the affine mapping introduced in 
Remark 15.71 The matrix Bt is defined in Remark 15.71 as well. And according to 
Remark [ST71 f(f) = f(f) for any / G A(T), where / is the image of / under the 
affine mapping from T to T. 

Definition 5.10. In the sequel, the phrase "for almost all" (parameters) will mean 
"for all except for a finite number" (of parameters). For example, a sequence x n in 
a topological space converges to x if, for every neighborhood of x, almost all values 
of x n belong to the neighborhood. 

Lemma 5.11. n.'~ T f : i? 1 (T;]R 2 ) — > 7 , f+1 A 1 (T) is a well-defined linear opera- 
tor for almost all t G [0, 1]. Moreover, 

diviJI.'~ f t & = n^diufcd) 

for any Cj G ^(TjR 2 ). 
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Proof. The linearity of the operator is obvious, for any isl. The point is to show 
that the operator is well-defined. For t = 0, the operator H.'" ^ Q reduces to the 

PB interpolation defined in [19] and, according to Lemma 9 and Lemma 20 from 
[T5] . is well-defined. Moreover, 

• fj, div x a}(x)i/j(x)diidi2 is a continuous functional of Cj, for any ip € V~,f\(T)/WL, 

• J^, w(x) T hi(x, t)dx\dx2 is a continuous functional of w, for any 1 < i < fcf 
and any igR, and 

• J-[Cj ■ n]?)ds is a continuous functional of Cj, for any ?) E "Pf(e)(e), and any 
ee A x (f). 

Therefore, in order to demonstrate that II ^ t is well-defined, it is sufficient to 
show that Cj = if Cj G V M A 1 (f) and II 1 " - Cj = 0. 

We take an arbitrary Cj G 'P f ^ 1 A 1 (T) such that ^--Tiffi = ^' According 
to the definition of T^A^T) and ((Ell), we know that cj € T^A^T). Set 
r = r(T). We denote by C(t,r) the matrix associated with conditions (15.41) and 
(|5.5[) . corresponding to a particular basis for 'P f ^ 1 A 1 (T') (the solution space), and 
a particular basis for V-^(T) /M.. We argue that matrix C(t,r) is non-singular for 
almost all t G [0, 1]. Notice that, for any r, det(C(t, r)) is a polynomial in £. Since 
U ?+i f o is well - derined > det(C(0, r)) ^ 0. So det(C(£, r)) is a non-zero polynomial. 
By the fundamental theorem of algebra, the polynomial has a finite number of real 
roots. We can conclude thus that l£<~ f t : H 1 {f;M 2 ) — > Vf +l A l (f) is well- 
defined for any t € [0, 1] except for the roots of det(C(i, r)). The number of roots 
is independent of the choice of the bases and depends only upon r = f(T). 

Since div*n^ f u E V f A 2 {f), for any Cj G H l {f; K 2 ), then IT 2 f div*n^ f Cj = 

div-tll 1 ' - . Cj. In order to show that div-o-D: 1 ' - Cj = Id 2 -div^a), it is sufficient to 

x f+l,T,t x f+l,T,t f.T x ' 

show that n 2 -div-v-a) = 0, for any Cj such that II.'" « Cj = 0. We have 

r,T x ' J f+l,T,t 

/ diVxWi/jdx = / div x w(x + c)dx = / div x wxdx + c / a) • fids = 0. 
if Jt Jdf 

Here IL'~ ^ t w = 0, ij> G "P?(f){T), c = ffipdji. The last equality holds due 

to the definition of n 1 ' - - and the fact that EL'~ ~ Cj — 0. We have thus 

f+l,T,t r+l,T,t 

divxll 1 " - Cj = n 2 ~div x w for any & G H^if; K 2 ). □ 

x r+l,T,t f,T x 

Theorem 5.12. Lei < G [0, 1] be any value for which Ilj"^ ^ t is well-defined. Then 

the operator T t : _ff 1 (T;IR 2 ) — 'Pr, 1 A 1 (T) zs well-defined as well, and we 

have the following result. 

There exist C > such thatjor T G Th, and ui G -ff^TjIR 2 ), the corresponding 
function Cj defined by 

(5.10) „(*(*)) = ^y-(x) 

belongs to H 1 ^ ;R 2 ) 7 and 

(5.H) n^- )T)tW( x(x)) = ^yntr.f/^) 
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Additionally, 

(5.12) divTi~~^ x Tt oj — II? T diuv 
and 

(5.13) ll n ?+i,T,t w IU 3 (T) < C||o;|| H i (T) 



Proof. It is easy to see that ui G iJ 1 (T;R 2 ). We define a linear isomorphism At 

from H 1 ^ ;R 2 ) to ii^TjR 2 ) by (A T ui)(x(x)) = , ?L M x). It is easy to see 

det(-Bx) 

that At is a linear isomorphism from Vf +1 h}(T) to 7 :, f+1 A 1 (T). 

By pulling back to T by (applied to both uj and H^Ti T t u) , we can see that 
(|5.7[) is the same as (|5.4p . (|5.8|l is the same as (|5.5p . and (I5.9[) is the same as (|5.6p . 
So we have (|5.11[) . As a consequence of (|5 . 11 1) , (|5 . 10[) and Lemma 15.111 we have 
(|5. 12|) . Inequality (j5.13[) is proved by using standard scaling techniques. □ 

Definition 5.13. We define a linear operator C t : [H^f)} 2 — ► P f+2 A°(f; R 2 ) by 
the following relations 



(5.14) J div i C t w(x)V'(x)dx = J div*w(x)^(x)dx V^e^iT)/! 

(5.15) / (C t w(x)) T h l (x,i)dx = / (w(x)) T h i (x,<)dx l<i<k? 
Jf Jf 

(5.16) A(C t t&) • n]f)ds = f[u ■ n]fjds Vfj G V m (e), Ve € A^f) 

(5.17) (\{C t u)-t]fids= f [u-i]fids V?) G P P( g)(e), Ve G Ai(f ) 

(5.18) C t w = at all vertices of T 

Here n, t denote the normal and tangent unit vectors along dT. 

Definition 5.14. For any T G Th, we define a linear operator Wx,t '■ H 1 ^; I 
V f+2 A- {T;R 2 ) by the following relations 



(5.19) / div(W T ,tW - u)(x)ip(x(x))dx = V^e V HT) (T)/R 
Jt 

(5.20) / (W T)t cj(x) - cj(x)) T S T T h. t (x(x), t)dx = 1 < i < k? 
Jt 

(5.21) [[(W T ,tW-u)-n]rids = Q Vry G V f{e) {e), Ve G A X (T) 



(5.22) / [(Wir,*w - w) • t]ijda = Q Vry G P P ( e) (e), Ve G Ai(T) 

J e 

(5.23) Wr.tw = at all vertices of T 

Here n, t denote the normal and tangent unit vectors along dT. As before, x = x(x) 
stands for the inverse of the affine mapping from the master to the physical element, 
defined along with matrix Bt in Remark 15.71 Also, according to Remark 15.71 
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r (/) = r(f) for any / 6 A(T), where / is the image of / under the affine mapping 
from T to T. 

Lemma 5.15. Operator C t : ff^fjR 2 ) — ► P ?+2 A (f ; R 2 ) is well-defined for 
almost all t G [0, 1] . 

Proof. Take an arbitrary f e [0,1]. Since w £ ^(TjR 2 ), then all right hand 
sides of (|5.14[) . (|5.15[) . (|5.16p . (|5.17j) are continuous functionals with respect to Cj. 
In order to show the well-dehnedness of Ct, it is sufficient to show that Cj = 
if Cj e 7 , f+2 A ('f;M 2 ), £> = at all vertices of T, and C t 6J = 0. According to 
conditions (15J6)) and (j5~T7| . E £> f+ 2A°(T; R 2 ). Define now C t : H x (f ;R 2 ) — > 
Pf=+2A°(T;R 2 ) by the relations 

(5.24) J div A C* t w(x)V3(x)dx = ^ div A w(x)^(x)dx V^(ic) e? f[f) (f)/i 

(5.25) / (C 4 w(x)) T h 4 (x,i)dx = / (w(x)) T h. ( (x, i)dx 1 < i < k? 
Jf Jf 

It is sufficient to show that operator C T ,t ■ H 1 ^^ 2 ) — > T f+2 A°(f ;R 2 ) is well- 
defined. 

Obviously, the right-hand side of conditions (15.241) and (|5.25l) are continuous 
functionals with respect to Cj. Set r = f(T) and denote by C(t, r) the matrix 
associated with the left-hand side of conditions (|5.24l) and (|5.25l) corresponding to 
some basis of Tf +2 hP(T\~R 2 ) (the solution space), and some basis of V F / f-, (T)/R. 

We need to show that the det (C(t, r)) ^ for almost all t £ [0, 1]. 

Since Cj vanishes on the boundary, we can integrate by parts (|5.24p without 
getting any boundary terms. So (|5.24[) is the same as 
(5.26) 

J^C t {Cj)(*)) T ■ grad^(x)dx = Jm±)) T • grad^(x)dx W- £ P f(T) (f)/R. 

Notice that ^ +2 A°(f;M 2 ) = x(x)[P?,f )_ X (T)] 2 where x is a third order poly- 
nomial vanishing along boundary dT and positive in the interior of T (a "bubble"). 
And notice that 

{hf,i(x, 1), • • • , hf, kf (x, 1)} © gradsVf (f > (T)/R 

={gr,i(*),-" ,g^(x)}©grad^ (f) (f) = [P f(f) _ 1 (f)] 2 . 

Consequently, det(C(l,r)) ^ 0. Since det((7(t, r)) is a polynomial in t, by the 
fundamental theorem of algebra argument again, det(C , (l,r)) vanishes at a hnitc 
number of roots only. 

We can conclude thus that operator C t : iJ^TjR 2 ) — ► 7 3 P+2 A (f ; R 2 ) is well- 
defined for almost all £ £ [0, 1]. Notice that the number of roots is independent of 
the choice of the bases, and depends upon f(T) only. □ 

Theorem 5.16. For any t £ [0,1] such that C t : iJ^fjR 2 ) — > P P+2 A°(f;R 2 ) 
is well-defined, there exists C > such that, for any T £ 7~h, operator Wx,t '■ 
H 1 ^;^ 2 ) — > P f+2 A°(T;R 2 ) is well-defined as well and, 

(5.27) \\curW T ,M\v<T) < CC^IMMt) + MhHt)) G H\T;R 2 ). 
Here hr is the diameter of T . 
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Proof. For any T £ Th, we define a linear isomorphism At from i? 1 (T;R 2 ) to 

i? 1 (T; R 2 ) by (Ara))(x(x)) = T . It is easy to see that At is a linear isomor- 

det(-fc>x) 

phism from P f+2 A a (f; R 2 ) to P F+2 A°(T; R 2 ). 

By pulling back to T by At (applied to both oj and Wr,tw), we can see that 
(|5.19|) is the same as (|5.14p . (|5.20|) is the same as (|5.15|) . (|5.21|) is the same as (|5.16p 
and (15.231) is the same as (|5.18l) . (|5.22l) becomes 
(5.28) 



[{A T l {W T ,tu)) 1 B T B T t]fids = / [Cj 1 B T B T t]fids V77 6 P f (e)(e), Ve e A X (T). 

^ e 

Notice that n and B^BTt are not parallel to each other. Since t T B^Brt ^ 0, 
t is not perpendicular to B^Bt^. Since n_Lt, then n is not parallel to B^Bt^. 
This means that conditions (|5.21l) and (|5.22l) imply conditions (j5 . 16[) and (|5.17p 
by the definition of At- So we can conclude that Wr,t is well-defined when Ct is 
well-defined, and 



(5.29) Wr,tw(x(i)) = A T C t u(&). 

Finally, inequality f|5 . 2T[) results from standard scaling techniques. 



□ 



Definition 5.17. Let r : A(T) — > Z+ be a locally variable order of discretization 
that satisfies the minimum rule. According to Lemma 15.111 and Lemma 15.151 there 
exists t~,fs e [0, 1] such that both II_'~ - and Ct - are well-defined. And the 

value of t-,fs depends on r(T) only. 

5.3. Projection operators on the whole mesh. According to Lemma 15.111 
Theorem [533 Lemma T5. 151 and Theorem 15. 16[ there exist {ti,t2> • " } C Z+ such 
that, for any T G Th, ^-f'+iTt- an d ^T,t HT) are well-defined linear operators. 



From now on, we rename II- 



l,- 



f+l,T,t HT) 



byn. 



-1,T' 



and Wrt-fxi by Wt- 



Definition 5.18. Wc define the following global interpolation operators, 

nl Th : L\n) —> v f A 2 (r h ), (n? Th u)\ T = n? T (u| T ) 
L 2 (r>;R 2 ) ^^,A 2 (r^;R 2 ), (np jTh U2 ) t )|t - (n 2 , T ( Ul | T ), n 2 , T ( U2 | T )) n 

-> P^AHTh), (TLl'- hTh u)\ T = TII^tHt) 

V f+1 A\T h ;R 2 ), (ni +liT ^)| T 



n 2 



u f+i,n 



ni+i 



: iJ^fijR 2 ) 



Til T12 
T21 T 22 



where 

Til 
T12 



- n f+i.r„ 



for all TeTfc. 



Oll\T 
0~12\T 
2 



T21 
T~22 



- n f+i,T„ 



0~2l\T 
C 2 2|t 



0"ll (T12 
C 2 i cr 22 



) — >• Pr +2 A°(T h ;R 2 ), (W Th u)\ T = W t (uj\t) 



Theorem 5.19. There exists C > such that 

rl,- 



l|n~ + i :rh er|| L 2 (s:2) < c\\o-\\ HHn) , \\nl: +1Th uj\\ L 2 {n) < c\\oj\\ H i {n) , 

for any a G # X (0;M) and uj G iJ^^R 2 ). Moreover, 



dMI~ +1 7- <7 = H 2 : j- h diva, 



div!i f ' +1Th uj 



n 



divui 
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Proof. This is an immediate result of Lemma 15.41 and Theorem 15.121 □ 

Definition 5.20. Let Rh denote the generalized Clement interpolant operator from 
Theorem 5.1 in [7], mapping i/^fijR 2 ) into V 1 h°{T h ; K 2 )- We define 

W h = W h (I-R h )+R h 

Theorem 5.21. There exists C > such that 

\\curW h 0j\\ L 2 {n) < C|M]ffi (n) G H\Sl;R 2 ) 

Operator Wh maps i? 1 (f2;R 2 ) into 7 , f=+2A°(T; IR 2 ) and satisfies the condition 

nif^w = nl>- 1Th W h u> e il 1 ^; M 2 ) 

Proof. We utilize Example 1 from [7J (in this case, Rh is the same as standard 
Clement operator) with uniform order equal 1 to construct operator Rh- Operator 
R h maps ff x (n;R 2 ) into ViA°(T; M 2 ) C TV^A^T; R 2 ). According to Theorem 5.1 
from [7J, there exists a constant c > such that, for any T € 7ft, 

(5.30) llw-i^lU^+HI^-^HIi^CT) < cMMIffHJrT) G H 1 ^; R 2 ) 

where -KTt stands for the patch of elements adjacent to T, .Kr = Ur'eTh T'nT^0 
As (7^)^ is regular, sup h sup TeTii #{T' E % ■ T' c ifr} < oo. We have 

|| curlWfcW |j L 2 ( n) <||curlW fc (7 - iih)w||i3(n) + HcurlRfcwH^n) 

^c^t^K 7 - R h)u\\ L *{n) + - Rh)u\\w(n) + \\cuT\R h L}\\ L 2 {n} ) 
<C||w|| fl i ( n) VwGff^R 2 ). 

The second inequality above holds by Theorem 15.161 and the third one by (|5.30[) . 

According to the definition of 11-^ r , and the definition of Wh, Hp'+i j-u = 
Ilb- hTh W h L>, for any uj G H^fyR 2 ). Notice that {I - R h )uo G ViA°(T h ) C H\Q). 
This implies 

s nh- ltTh uj = n.h- ltTh w h u}. a 

6. Stability of the finite element discretization 
We need the following well-known result from partial differential equations, see 

e.g. m- 

Lemma 6.1. Let be a bounded domain in R 2 with a Lipschitz boundary. Then, 
for all fi G L 2 (fl), there exists n G i/ 1 (fi;]R 2 ) satisfying diuq = /i. //, in addition, 
J Q H = 0, then we can choose n G i? 1 (fi). 

The main result of this paper for afiine meshes is Theorem 16.31 below. Its proof 
follows the lines of Theorem 9.1 from [14], Theorem 7.1 from [6] and Theorem 11.4 
from [5]. The main difference is the use of operator Wh in place of the operator 

nr 2 m m- 
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Lemma 6.2. There exists c> such that, for any (tu , y) G Vf A 2 (Th) xVf A 2 (%; R 2 ), 
there exists a G Vf+iA 1 (7h',^ 2 ) such that 

diva = /i, —II 2 j- h S\<J = ui 

and 

\W\\ H (div,n) < c{\\n\\ L 2 {n} + ||w|| L 2 (n) ) 
Here, the constant c depends on swp h sup Tg q- f(T). 

Proof. (1) By Lemma lOl we can find r\ G H 1 ^; M) with div?? = /i and H^Hi^si) < 
c||/i||z,3(n)- 

(2) Since w + IT? h SiUl +1 h n G L 2 (£l), we can apply Lemma I5TT1 again to find 
r G iJ^^R 2 ) with divr = uj + II? ^Sifti^^ such that 

IM| ff i(Q) < c(||w|| L 2 (n) + ||II? h Sifl} +lth r]\\ L 2 (n) ) 

(3) Define ct = cwclW h T + ftl +1/t ?7 £ Vf+iA^T^M. 2 ). a G P f+ iA 1 (%] R 2 ) by 
Lemma 13.51 

(4) From Step (3), Theorem EH Step (1), and the fact that ft? ^rf? ^ = ft? jh , 
we obtain 

diver = divnl +1A ?7 = ft-^div?? = ft? h n = /i. 

(5) Notice that divw = — S'lCurlw, for any w G i? 1 (fi;R 2 ). Here operator S\ is 
defined in ([i~10]) . Then II? ^diviri^ h cj = -II? ^Sicurlu;, for any u G i^^R 2 ). 

(6) Also from Step (3), -II? h S 1( T = -II? curler - n? h SiiF= +1 ^ Apply- 
ing, in order, Step (5), Theorem I5.2H Theorem 15.191 Step (2), and the fact that 
n^n? h = H? h , we have 

-n? ih 5icuriWhT = n},, h dwiil'~ lh w h T = n? ^divn^^r = n 2 Vl divr 

= n L( w + U f,h^M+i,hV) = u + U f.hSifll +lh r,. 

Combining, we have —II? h S\a = uj. 

(7) Finally, we prove the norm bound. From the boundedness of Si in L 2 , 
Theorem 15. 191 and Step (1), we obtain 

||n? fc 5'ini +1)fe 77|| z ,2 (n) < c||s , iiii +1 ^77|| i 2 (n) < c||ni +1)fe 77|| i2(n) 

^ c \\v\\ H ^Q) ^ c llMlli2(o)- 

Combining the result with the bound in Step (2), we get ||r||fl-i(m < c ( II ^ II z> a fr2) + 
HmIIl 2 ^))- Theorem 15.211 then yields 

llcurlWfcrllz^) < c||r|| H i ( n) < c(|H| L2(n) + |HIl 2 (q))- 

From Theorem 1 5 . 1 91 and the bound in Step (1), 

W^lr+i.hVh^n) < ch|| ff i(o) < c\\ f i\\ L 2 {n) . 

In view of the definition of a, these two last bounds imply that ||cr||L 2 (f2) < c (ll £ *'||i, 2 (f2) + 
IHU 2 (n)) ; while ||divo-|| L 2 (a ) = \\[i\\ L 2(n) by Step (4). This finishes the proof. □ 
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Theorem 6.3. There exists c > such that, for solution (o~,u,p) of elasticity 
system il.8\) . and corresponding solution (o~h,Uh,Ph) of discrete system fi4-l\ l, we 
have 

Ik - VhWiiidiv.n) + \\u - u h \\ L 2( n) + \\p -ph\\L 2 (n) 

<cinf[||(7 - r\\ H ( diVi n) + II" - w|U a (fi) + lb _ 9lU 2 (n)L 

where the infimum is taken over all r G Vf+iA 1 (Th, R 2 ), v G VfA 2 (Th,M. 2 ), and 
q eP f A 2 (T h ). 

Proof. We need to show that conditions (|4.3[) and (|4.4I) are satisfied. Condition (I4.3|) 
follows from the fact that, by construction, dxvVf+iA 1 (%, M 2 ) C V? A 2 (%] R 2 ), and 
the fact that A is coercive. Condition (|4.4j) comes from Lemma T6.2I This finishes 
the proof. □ 

7. Curvilinear Meshes 

In practice, meshes generated by CAD software are usually curvilinear. In the 
following sections, we generalize our result for affine meshes - Theorem [673] to curvi- 
linear meshes in the sense of asymptotic ft-stability. 

7.1. Mesh regularity assumptions. 

Definition 7.1. (Curved triangle) A closed set T C R 2 is a curved triangle if there 
exists a C^-diffeomorphism Gt from reference triangle T onto T . This means that 
G T is a bijection from f to T such that G T G C 1 (f) and G^ 1 G C 1 ^). We assume 
additionally that det (DGt(x)) > for any x G T. 

We represent Gt in the form 

(7.1) G T = G T + $t, 

where Gt : A — ^ i?T& + br, P>t = DGt(P) with p being the centroid of T, and $t 
a G 1 -mapping from T into M 2 . The images of edges and vertices of T by Gt are 
edges and vertices of T respectively. We denote A,(T) = Gr(Aj(T)),i = 0,1,2, 
and A(T) = G T (A(f )). 

Definition 1 7. II is practically identical with the definition of curved finite elements 
introduced in [?]■ 

Definition 7.2. A curved triang le T is of class G fe , fc > 1, if the mapping Gt G 
C k (T). Similarly, a curved triangle T is of class C ' ,k > 1, if the mapping Gt G 
G fe ' 1 (T). 

We define 7^ to be a finite set of curved triangles T, where h denotes the maximal 
distance between two vertices of T G Th ■ We define vertices of Th to be vertices of 
T G Th, and we define curves of Th to be edges of T G Th- We assume that any 
edge of T G Th is either an edge of another curved triangle in Th, or part of the 
boundary of Th- 

Each curve of Th is parametrized with a map from the reference unit interval 
into R 2 , 

[0,1] 3S4 x e (s) G M 2 
The parametrization determines the orientation of the curve. 
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Let £(s) be the local parametrization for a particular edge of a curved triangle 
T £ Th, occupied by a curve e of Th- This means that £(s) is an afhne mapping 
from the reference interval onto an edge of the reference triangle T, whose image 
under the mapping Gt is exactly the particular edge of T. We can choose £(s) so 
that Gt(C(s)) has the same orientation as x e (s). 

Definition 7.3. (C°-compatible mesh) Th is called a C°-compatible mesh if, for 
any curve e and any curved triangle T which contains e as an edge, there is a local 
parametrization £(s) of e satisfying 

G T (C(s))=x e ( S ). 
The concept is illustrated in Fig [TJ 




Figure 1. Compatibility of edge and triangle parametrizations. 
We denote 

(7.2) c h := sup((sup|| J D$T(x)||)||i? T 1 ||) 

Ten ief 

For each T £ Th, we define T = Gt (T) . We denote by hx the diameter of T and 
by pt the diameter of the sphere inscribed in T. 

We define A z {T h ) = [j TeTh A,T, and A(%) = U TeTh AT. 

Remark 7.4. In order to simplify analysis, compared with [7], our definition of (j7.2[) 
replaces \\D$ T ■ B^\\ with the upper bound H-D^tII • H-B^ 1 !!- 

Definition 7.5. The family {Th)h of C°-compatible meshes is said to be regular if 

sup sup hr/pr = & < oo, and lim Ch = 

h Ten 

We show the construction of (Th)h of C°-compatible meshes in Appendix [XJ 

Lemma 7.6. There exist C\,ci > such that, for any triangle T, 

ci\\Bt\\ ■ WB^W < h T /p T < c 2 \\B T \\ ■ WB^W, 

where x = Bt~X- + br is the affine homeomorphism from T to T. 

Proof. Ht I Pt < C2||-&r|| • Hi?^ 1 !! comes from the geometric meaning of singular 
values of matrix Bt- ci\\Bt\\ ■ H-B^ 1 !! — ^t/pt is a consequence of Theorem 3.1.3 
in [9]. □ 
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Lemma 7.7. Let family (Th)h be regular. Then, for any indices i,j,k,l € {1, 2, 3}, 
we have 

lim sup sup U 7^-, r = 

h^o TeTh±ef det(B T ) 

Proof. For any T e Th and any x G T, we have 

Since H-B^ll ' II b t|| > 1, _ x < ||S T ||. Consequently, 

ll-^T II 

' W toW W)M | " (l|2? * r(i)l1 ' llVlDHSrllVldetWI. 
Since sup^ supyg-j-^ hx/pr — cr < oo, ||Bt|| 2 /| det(Sr)| < or 2 with c > 0. 

Since hm^oCf, = 0, we have hm ft _>o sup TeTh sup^ | ) 1 = °- 

□ 

Lemma 7.8. If a family (7~h)h is regular, then we have 

lim sup sup WB^DGrix))- 1 - I\\ = lim sup sup UpGrCx))" 1 ^ - I\\ = 0. 
^reTiigf Ten xef 

Proof. For any T e Th and any x S T, we have 

H^DG^x))" 1 -!]] = HBt^t + ^tCx))- 1 -/]! = iki + ^t^b^ 1 )- 1 -/]]. 

Since — ► as h — » 0, we have lim7 l _ J ,o su PTeTh su P x gf II-St(-DGt(x)) _1 — I\\ = 0. 
Proof of the second property is fully analogous. □ 

Lemma 7.9. If a family (Th)h * s regular, then we have 

lim sup sup\det(DG T (x)B^)-l\ = lim sup sup | det^T^Grlx))" 1 )-!! = 0. 

Proof. Since lim^o Ch — 0, lim^o su PTeT h su Pief I det(JJGT(x)i?y 1 ) — 1| = 0. 
By Lemma ITTS1 lim/j_yo su PTeTfc su Pief I det(BT(-DGT(x)) -1 ) — 1| = 0. □ 

8. Finite element spaces on curvilinear meshes 

We begin by introducing the relevant finite element spaces on any curved triangle 
T by the pull back mappings associated with the inverse of Gy, where Gt maps 
from the reference triangle T to T, see Definition 17.11 Then, we define the finite 
element spaces on a whole mesh Th by "gluing" the finite element spaces on curved 
triangles. 

8.1. Finite element spaces on a curved triangle. Let T be a curved triangle 
from Definition 17.11 with Gt denoting the corresponding G^diffeomorphism from 
T to T, x = Gt(x). We begin by introducing formally the mapping f from A(T) 
to Z_|_ specifying the local order of discretization. 

Definition 8.1. We denote by f a mapping from A(T) to Z + such that if e, / 6 
A(T) and e C / then f(e) < r(f). With the same symbol f we denote the cor- 
responding mapping from A(T) to Z+, r(f) ■= r(f) for any / G A(T), where 
f = G T (f). 
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We define now the following FE spaces on T: 
Definition 8.2. 

T f A°(T) := {u(x) : u(x) G T f A°(f) where u(x) = u(x)}, 

P f A x (T) := Mx) : w(x) G T f A\f) where w (x) = 1— - J)G T (x)a;(x)}, 
Vf h}(T) := { W (x) : £(x) G PrA^f ) where w (x) = — ^1— - £)Gfr(x)fl(x)}, 

VfA 2 (T) := {«(x) : fi(x) € 7> f A 2 (T) where u(x) = ^— 1-— 6( x )}, 

7> ? A°(T;1R 2 ) :={(ui,u 2 ) : Ul ,u 2 GP f A°(T)}, 
7^A 2 (T;1R 2 ) := {(ui.ua) : ui,u a G PfA 2 (T)}, 



»2x 



Oil (Ti2 
021 022 



(ffii,<T 12 ) T ,( ( T 21 ,a 22 ) T G A 1 (T) 



Remark 8.3. Since Gt : T — > T is a C 1 -diffeomorphism with det(-DG;r(x)) 7^ 0, for 
any x G T, the formulae in Definition 18.21 arc well-defined. The mappings used in 
the definition are the standard pull back mappings for differential forms A , A 1 , A 2 , 
see e.g. formulas (2.24), (2.26), (2.27) in [15] . 

Lemma 8.4. For any edge e G Ai(T), Zei £(s) ^ e ^ e local parametrization for e 
discussed in Definition \7.3\ i.e. the affine mapping from the reference interval onto 
e G Ai(T). We have then, 

V f A°{T)\ e = {u(x) w/iere x G e : u(G T (C(*))) G Pf(fi)(e)}, 
TVA^T) • n| e = {u(x) w/iere x G e : u(G T (((s)))\\D(G T o ()(a)\\ G V m (i)}, 
V~A\T) • n| e = {«(x) where x G e : W (G T (C(s)))P(G T o C)(s)|| G P f(g) _ 1 (e)}. 

In addition, the above equalities do not depend on the choice of the orientation of 
the local parametrization £(s). 

Proof. These are trivial observations on pull back mappings and their restrictions 
to edges. □ 

8.2. Finite element spaces on a G°-compatible mesh. Let 7ft, be a C°-compatible 
mesh from Definition 17.31 We extend the f mapping to a global map defined on 
A (7ft) with values in Z + such that if e C /, then f(e) < f(f). 

Definition 8.5. We put flh ■= UtsT, ^ 

GA°(7ft) := {u G i/ 1 (ilft) : u is piece-wise smooth with respect to 7ft}, 
GA 1 (7ft) := G Z7(div, fi^) : u> is piece-wise smooth with respect to 7ft}, 
GA 2 (7ft) := {u G L 2 (Qh) : u is piece-wise smooth with respect to 7ft}. 
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We define 



V f A°(T h ) := {u G CA°(T h ) : u\ T G P P A°(T),VT € TJ, 
PfA^^) := G CA^Tfe) : u\ T G V f A\T),\/T G 7"h}, 
VfA 1 (%) := {lo G CA 1 (7h) : w|t G V? A l {T),\fT G 7^}, 
P f A 2 (T ft ) := G CA 2 (T h ) : u|r 6 P=A 2 (T),VT G T h }, 
P f A°(r,;R 2 ) := [T 3 ? A°(7/i)] 2 , 'PfA 2 (Th; R 2 ) := [7> f A 2 (7;)] 2 , 



Remark 8.6. According to Lemma [8.41 and the fact that 7ft is C°-compatible, we 
can conclude that V f A°{T h )\ T = Pf A°(T), V? A 1 (%)\t = VfA 1 (T),VfA 1 (%)\T = 
^A 1 (T),7'fA 2 (rh)|T = VfA 2 (T), for any T G 7*. For standard (not curved) 
triangulations Th, spaces V? A°(Th), V? A 1 (7/ 1 ), 7 3 f TA 1 (7ji 1 ), V? A 2 {Th) coincide with 
those analyzed in [T9] . 

Lemma 8.7. 



9. Preliminaries for the proof of stability for curvilinear meshes 

We begin by recalling our assumptions on the domain and meshes: Q is a (curvi- 
linear) polygon and it is meshed with a family (7~h)h of C°-compatible meshes of 
class C 1 ' 1 . For any mesh Th, mapping f : A(Th) — ► ^+ defines a locally variable 
order of discretization that satisfies the minimum rule. The maximum order is 
limited, i.e. sup^ sup T67 - h f(T) < oo. 

As in Section [5j we shall design operators LI 2 T , II- j. and Wt in order to 
make the left and the right diagrams in (|4.5[) commute. In the case of affine meshes 
(Section [5]), we define LI-^ T and Wt on physical triangles, then "pull them back" 
to the reference triangle. We designed Hf'+i T and Wt in such a way that their 
"pull backs" are Yll'~ ~ and Ct_ . respectively ( Definition 15 . 1 T|) . In the 

case of curvilinear meshes, the approach to Lip T and II-C^ T is similar to that in 
Section [5j But the treatment to Wt is quite different. In order to make the right 
diagram in (|4.5p commute (we need the commutativity on the physical meshes), 
we incorporate the construction of II -i^ T into the definition of Wt- Then we use 
a special mapping such that the "pull back" of Wt is a projection operator to 
Vf+2A°(T; R 2 ). In general, this kind of "pull back" of Wt to the reference triangle 
does not coincide with Ct f{T) from Definition 15.131 This means that, in general, 
Wt is not necessarily well-defined. But we manage to prove that the "pull back" 
of Wt will converge to C tf(T) as h — ¥ 0. So Wt becomes well-defined when h 
is small enough. The work presented here illuminates the substantial difference 
between the stability analysis for curvilinear meshes and that for affine meshes. 
This is in particular reflected in the proof of Theorem I9.15[ which uses heavily 
Lemmas 17.717.817.91 





PfA^Th) c V^ +1 A\T h ) c Vr+xA l {T h ), 
divVf+iA^Th) C T f A 2 (T h ), curtP f+1 A Q (T h ) C VfA 1 



(Th). 



Proof. The proof is straightforward. 



□ 
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(9.1) / (Il£ x u(x) - n(x))^>(x(x)) dx = Vip £ V f{T) {T) 



In order to make the paper more readable, most proofs of results discussed in 
this section have been moved into Appendix [Cj 

Definition 9.1. For any T G Th, we define a linear operator II? T : L 2 (T) — > 
V f A 2 (T) by the relations 

S 2 

it 

Above, x(x) signifies the inverse of the element map x = Gt(x). 

Definition 9.2. Operator II 2 - : L 2 (f) — > V f A 2 {f) will denote the L 2 -projection 
in the reference space, 

(9.2) /(ie t fi(*)-fi(x))^(*)dft = V4>GVr iT )(f) 

JT 

Remark 9.3. Operator II? T is a weighted L 2 -projection in the physical space. For 
a regular triangle (affine element map), the jacobian is constant, and 11? T reduces 
to the standard L 2 -projection in the physical space. 

Lemma 9.4. For any T € Th, and arbitrary u(x) £ L 2 {T), we define u(x.) by the 
relation: 

U(X(X)) = dct(DG T (x)) 

TTien w(x) G L 2 (T), and 

n 2 f «(x) 

^^ (X(X)) = detp GT (x))) 
Above, x(x) signifies the element map x = Gr(x). 



Proof. Proof follows immediately from Lemma IB.6I and the definitions of the two 
projections. □ 

Lemma 9.5. For any e > 0, there exists S > such that, for any h < 5 and 

Te%, 

||n? ]T M - Pf,Tu\\ L 2(T) < e\\u\\ L 2( T ),\/u G L 2 (T). 

Here Pf t T is the standard L 2 -projection onto Vfk 2 (T). 

Proof. Please see Appendix [Cj □ 
9.1. Projection Based Interpolation onto 7 , f+iA 1 (7/ t ). 

Definition 9.6. For any T G Th, we define a linear operator Ili +1 T : 7J 1 (T; R 2 ) — !► 
Vf+iA 1 ^) by the relations 

(9.3) / div(ni +1T w-w)(x)^(x(x))dx = V$ G V« T) (f )/K 
Jt 

(9.4) / (nl +1 tw(x) - w(x)) T J DG T (x)" T curlx^(x(x))dx = V<£> G Pf(T)+ 2 (T) 

JT 

(9.5) 

/ [(11^! T w-w)(x e (s)).n(x e (fl))]i)(«)||x e (s)||ds = Vt) G 7> f(e)+ i([0, l])Ve G A X (T) 
i[o,i] 
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Here x = Gt(x) for any x G T, x e (s) : [0, 1] — > e is the parametrization of e, and 
n is a unit normal vector along e (the choice of its direction does not matter) . 

Definition 9.7. (Projection Based Interpolation operator onto Vf+\h}(T)) We 
define a linear operator II;. ^ : H(T) — > Vf+iK l (T) by the relations 



r(T) 



(9.6) jf div*(n^ 1)# w-w)^f* = vi>ev H 

(9.7) y^ +lf w-w)-curl^dx = V^eP f(f)+2 (f) 

(9.8) f(nl +lf ij-uj)-hf)ds = \ffjeV m+1 (e)yeeA 1 (f) 

J e 

Remark 9.8. The operator n* +1 ^, is the Projection Based Interpolation operator 
onto V f+ iA 1 (T) defined in [T5]. The operator n~ +1 T is defined by the pull back 
mapping from T to T. 

Lemma 9.9. For any T G 7h, (Miy w G [H l (T)] 2 , we define Cj by the relation 



DGt(x) , m 
" (X(X)) = d«xj) uW 



TTien w(x) G H{f), and 



Proof. Please see Appendix [Cj □ 

Lemma 9.10. For any T G Th, and any uj G [i? 1 (T)] 2 J we Ziawe Hf T divuj = 
divn.f +1T u). 

Proof. Please see Appendix [Cj □ 
Lemma 9.11. There exists 6 > and C > suc/i t/iat, /or any h < 8, we have 

||n^ T w|| L 2 (T) < C||w|| ff i ( r) VT G 7h, ^ £ H X ^T; R 2 ) 
-For affine meshes, the inequality above holds for any h > 0. 

Proof. Please see Appendix [UJ □ 

9.2. Modified Projection Based Interpolation onto Vf +1 A (T) and modi- 
fied operator W onto P f+2 A (T). 

Definition 9.12. (PB interpolation operator onto Vf +1 h}(T)) For t~,f^ given 
in Definition 15.171 and for any T G 7fc, we define a linear operator 11-^ T : 
H^T^ 2 ) — > Pp +1 A 1 (T) by the relations 

(9.9) / div(np\7 1T w-a;)(x)^(x(x))dx = G P f (T) (f)/K 

(9.10) / (nlf 1)T w(x) - w(x)) T ^G T (x(x))- T n l (x(x), t f(#) )dx = 1 < * < 
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(9.11) 

f [(III'-, T o;-a;)(x e ( S ))-n(x e ( S ))]r}( S )||x e ( S )||d S = Vr) G P f(e) ([0, 1]), Ve £ Ai(T) 

^ [0,1] 

In the above, x = x(x) signifies the inverse of the element map. 

Theorem 9.13. The operator Ul'^ T : ii^T;^ 2 ) — > Vf +1 A l (T) is well-defined, 
and we have the following result. 

divnl'+ 1)T u) = IJ$ tT divu Vw eff^TjK 2 ). 
There exist 5 > and C > such that, for h < 5, T G Tu, and ui G H (T; M 2 ), 

ll n f+l,T w IU 2 (r) < C\\uj\\ h i {t) 

Proof. The proof is analogous to that of Lemma [9791 Lemma r9.10[ and Lemma[9Jj] 

□ 

Definition 9.14. For t f ^ given in Definition 15. 171 and for any T G 7ft, we define 
a linear operator W T ■ F^TjM 2 ) — > 7^+2 A (T; R 2 ) by the following relations 



(9.12) y div(Wrw - uj)(x)tp(x(x))dx = VV> € Pf(r)(T)/R 

(9.13) y (W T w(x)-a;(x)) T L)GT(x(x))- T n i (x(x),^ (f) )fix = 1 < i < k f 
(9.14) 

[(W r r,tw-w)(x e («))-n(x e ( a ))]^(a)||x e («)||dfl = Vr) G P f(e) ([0, l]),Ve £ A X (T) 



[0,1] 

(9.15) 

[(Wr,tw-w)(x e (a))-t(x e (s))]^(«)||x e («)||dj = V?) e Pf (e )([0, l]),Ve G Ai(T) 

[0,1] 

(9.16) W T ,tUJ = at all vertices of T 

Here n, t denote the normal and tangent unit vectors along dT. 

Theorem 9.15. There exist 5 > and C > sucft i/iai, for any h < S, T G 7ft, 
operator W T ■ H l (T;R 2 ) — ► P F+2 A°(T; M 2 ) is well-defined and, 

\\curW T u\\ L i {T) < C(h^\\ij\\ L 2 {T) + ||w|| ff i (T) ) G i/ 1 (T; M 2 ) 

for affine meshes, Wt is well-defined and the above inequality holds for any h > 0. 

Proof. Please see Appendix O □ 

9.3. Projection operators on the whole curvilinear meshes. 

Definition 9.16. Similar to Definition 15.181 We define the following global inter- 
polation operators, 

n 2 ^ : L 2 (n) —> VrA 2 (T h ), (uj iTh u)\ T = n? T (u| T ) 
iil Th -.L 2 (n-R 2 ) -^v f A 2 (T h ;R 2 ), (n? T( >i, U2 ) T )| T = (n 2 , r ( Ul | r ), n 2 , T (u 2 | T ))" 
n£i i7i : h^k 2 ) — > v^a^Th), (n£ ltTh u)\ T = iferMr) 
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where 

Til 
T12 



- n f+l,Th 



for all Te%. 



C12|t 

2\ 



T21 
T"22 



= m 



V+l,^ 

Pf+ 2 A (Th;]R 2 ), 



0"2l|T 
02 2 |t 



0"ll 0"12 
0"21 ^22 



(Wr h u)\T = W t (oj\t) 



Remark 9.17. Since (Th)h is C°-compatible, operators II f ^ 1 Th , fl} +1Th and Wr h 
are well-defined. 

Theorem 9.18. For any e > 0, i/iere exists S > smc/i i/ia^ /or any < 8, 



in 



^,r h u||L2(a) < e||«|U 2 (n) Vu e L 2 (fi) 



-Here Pf.T h * s ^ e standard L 2 -projection onto Vfh 2 {Th)- 

Proof. This is an immediate result of Lemma 19.51 □ 

Theorem 9.19. There exist 5 > and C > suc/i i/iai, for any h < 5, we have 

ll n f + i,rh cr lli 2 (n) < C\\a\\ H i (n)l Hn-^ Th a;|| L 2 (0) < C\\oj\\ H i {n) , 

for any a £ H^f^M) andu G ii/^iljR 2 ). For affine meshes, the inequalities above 
hold for any h > 0. Moreover, 



divfl} +ljTh ( 



Il| j- diva, 



divHl' +ltTh ui 



II? t- divu) 



Proof. This is an immediate result of Lemma 19.111 and Theorem 19.131 



□ 



Definition 9.20. Let Rh denote the generalized Clement interpolant operator from 
Theorem 5.1 in [7], mapping ff^R 2 ) into VxA (T h ; M 2 ). We define 

W h = W h {I - R h ) + R h 

Theorem 9.21. There exist S > and C > such that, for any h < 5, 

\\curW h uj\\ L 2 in) < C\\u\\ mm e H'iQ-M 2 ) 

For affine meshes, the inequality holds for any h > 0. Operator Wh maps H 1 (Q; R 2 ) 
into 'P,~+2A°(T; R 2 ) and satisfies the condition 

n&i,r h w = ^l+i,n w '^ Vcj e ffl ( n ; r2 ) 

Proof. We utilize Example 2 from [7] with uniform order equal 1 to construct 
operator Rh- Since (7~h)h is C -compatible, and Ch — >■ as h — > 0, operator i?^ 
maps ff 1 ^;^ 2 ) into 7\A°(T;R 2 ) C V? + 2A°(T; M 2 ). The the proof will be the 
same as that of Theorem 15.211 □ 

10. Asymptotic stability of the finite element discretization on 

curvilinear meshes 

Lemma 10.1. There exist 6 > and c > such that, for any h < S and any 

6 V f A 2 {%) x V f A 2 (T h ;M. 2 ), there exists a € Vr+iA 1 (%;M 2 ) such that 



diva = /!, — ri~7-,Si<7 = 



and 



\W\\ H (div.n) < c(]HU=(fi) + IM|z,3(n)) 
flere, t/ie constant c depends on sup h su l?TeTh ^(-^)- ^ 7 ' or a ffi ne meshes, the inequal- 
ity above holds for any h > 0. 



MIXED FEM FOR 2D LINEAR ELASTICITY 



27 



Proof. The proof is the same as that of Lemma 16.21 □ 

Theorem 10.2. There exist S > and c > such that, for solution (a, u,p) of elas- 
ticity system il.8\) , and corresponding solution (ah,Uh,Ph) of discrete system \4-ty > 
we have 

\W ~ o-h\\H(div,n) + \W - u h \\ L 2 {n) + \\p-ph\\mn) 
<einf[||cr - T\\ H{div ^) + \\u - u||i=(n) + lb - g|U=(Q)], 

where the infimum is taken over all r 6 Vf+\A 1 {7h, K 2 ), v S V? A 2 (Th, R 2 ), and 
q £ V?A 2 (Th), for h < S. For affine meshes, the inequality holds for any h > 0. 

Proof. We need to show that conditions (|4.3|) and (|4.4[) are satisfied asymptotically 
in h. Condition (14.31) follows from the fact that, by construction, divVf+i A 1 (Th, R 2 ) C 
P?A 2 (Th', R 2 ), and the fact that A is coercive. 

We turn now to condition (14.41) . According to Lemma 110.11 there exist 5 > 
and c > such that, for h < S and G T f A 2 (T h ) x T^A 2 ^; M 2 ), there 

exists (7 e Vf+\A l {Th] K 2 ) such that diver = /i, —Hp j- h S\o- — u, and ||0"||H(div,f2) < 
c(IHU*(n) + ||w|| L2(n) ). We have then 

(divcr,^) - (Si<t,u) =(divcr,/i) - (n 2 /7 - h S'icr,a;} + ((LT 2 ^ - P ftTh )Sia,u>) 

>c\W\\H(divM)(\\^\\L^n) + N|z,3(n)) + ((n 2 , Th -Pf,rJSio-,w). 
According to Theorem 19. 181 for sufficiently small h, 

\{Q%,T h -Pf,T h )Sia,u)\ < ~\\o-\\mn)ML*(n) 
So, asymptotically in h, we have 

(divcr,^) - (Sia,u) > -||o-|| ff ( div ,f})(IHU 2 (n) + 11^ )- 

For affine meshes, H 2 j- reduces to the standard L 2 -projection. The inequality 
above holds then for any h > 0. This finishes the proof. □ 

11. Numerical experiments 

We continue numerical experiments initiated in [19] where we investigated rates 
of convergence for uniform /i-refinements in presence of non-uniform polynomial 
order p, and tried out p-adaptivity (with no underlying theory at presence) . The ex- 
periments confirmed the optimal ^.-convergence rates and indicated the p-convergence 
as well. 

Following [T5] , we consider the L-shape domain, and use a manufactured solution 
corresponding to the exact solution of the homogeneous equilibrium equations (zero 
volume forces), and the corresponding unbounded L-shape domain extending to 
infinity. The manufactured solution is designed in such a way that both p and all 
stress components have the same singularity at the origin characterized by term 
r -o. 39596 w h ere r j s the distance to the origin. 

Experiments presented below focus on ft,-adaptivity for meshes with uniform 
order p, and the ultimate goal of this research - the /ip-adaptivity. We investigate 
both affine and curvilinear meshes. We use the standard "greedy algorithm" for 
/i-rcfincmcnts, and the two-grid /ip-algorithm for fop-refinements, see [12] for details. 
The /ip-algorithm is based on the standard Projection Based (PB) interpolation. 
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All convergence plots are displayed on log-log scale, error vs. number of degrees- 
of-freedom (d.o.f.). The error is always measured in terms of percent of the total 
norm of the solution. 

Uniform /i-rennements on affine meshes. We begin with a verification of sta- 
bility on uniform meshes. Fig. [2] presents the L-shape domain with an initial mesh 
of six elements. Fig. [3] displays the actual approximation error compared with the 
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Figure 2. The L shape domain with initial mesh. 

best approximation error for a sequence of uniformly refined meshes of zero ordeiQ. 
As expected, the two curves are parallel to each other. We repeat now the same 




10° ' 1 — ■ — 1 — ■ — I 

10 10 10 10 10 

number of d.o.f 

Figure 3. Uniform /i-refinements for p — 0. FE error vs. best 
approximation error. 

experiment starting with an initial mesh of elements with order varying from zero 
to four, shown in Fig. |4] The corresponding FE error is compared again with the 
best approximation error in Fig. [5j The two lines are again parallel to each other 
with the slope determined by the lowest order elements in the mesh (same as in the 
first example). 



'We always refer to the order of approximation for the displacement. 
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Figure 4. The L shape domain with initial mesh. 




number of d.o.f 



Figure 5. Uniform h- refinements for a mesh of variable order, 
error vs. best approximation error. 



FE 



Adaptive h- and ^-refinements on affine meshes. We continue now the 
numerical verification of stability for non-uniform meshes resulting from h- and 
/ip-refinements. Unfortunately, we face a slight discrepancy between the presented 
theory and the numerical experiments as the code is using 1-irregular meshes which 
we have not accounted for in our theoretical analysis. 

Fig. |5] presents convergence history for adaptive /i-refinements for meshes of 
order p — 0, 1, 2 and adaptive /ip-refinements starting with the mesh of zero order. 
The approximation error is compared with the best approximation error computed 
on the same meshes (generated by the adaptive algorithm). Results for /i-refined 
meshes of order p = 1,2, and the /ip-refinements confirm the stability result. The 
result for the lowest order elements, though, reflects some loss of stability. The 
only possible explanation that we have at the moment, is the use of meshes with 
hanging nodes. 

Finally, in Fig. we compare the convergence history for all tested refinements. 
The /ip-adaptivity produces the best results although, in the presented range, the 
rate seems to be still only algebraic. 

Adaptive h- and /ip-refinements on curvilinear meshes. We use the same 
manufactured solution on the circular L-shape domain shown in Fig. [8j The six 
triangles in the initial mesh are parametrized using the transfinite interpolation 
technique, see [12] , p. 201, for details. The FE error is compared again with 
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(a) Adaptive /i-refincments, p = 0. (b) Adaptive ft-refinements, p = 1. 




(c) Adaptive /i-refincmcnts, p = 2. (d) Adaptive /tp-rcfinements. 

Figure 6. Comparison of FE error for adaptively refined affinc 
meshes with the best approximation error. 



the best approximation error for a sequence of h- adaptive meshes and p = 0,1,2, 
in Fig. HI The results indicate again a slight loss of stability for the elements of 
lowest order. Finally, in Fig. I10| we compare the convergence history for all tested 
refinements. The /ip-adaptivity delivers again the best results but the exponential 
convergence is questionable. This may indicate that the use of standard Projection- 
Based interpolation operator in the reference domain is not optimal. 



12. Conclusions 

We have presented a complete /i-stability analysis for a generalization of Arnold- 
Falk-Winther elements to curvilinear meshes of variable order in two space di- 
mensions. The stability analysis for both generalizations: variable order elements, 
and curvilinear elements proved to be rather non-trivial. The case of variable el- 
ements has been tackled with a novel logical construction showing the existence 
of necessary interpolation operators rather than constructing them explicitly. The 
presented construction departs from operators used by Arnold, Falk and Winthcr 
and modifies the Projection-Based (PB) interpolation operators as well. 
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Figure 8. A circular L-shape domain with an initial mesh of six elements. 

The analysis of curvilinear meshes for elasticity differs considerably from that 
for problems involving only grad-curl-div operators. Piola maps transform gradi- 
ents, curls and divergence in the physical domain into the corresponding gradients, 
curls and divergence in the reference domain. Consequently, problems involving the 
grad, curl or div operators only (e.g. Maxwell equations or the mixed formulation 
for a scalar elliptic problem) can be reformulated in the parametric domain at the 
expense of introducing material anisotropics reflecting the geometric parametriza- 
tions. This is not the case for elasticity where the strain tensor (symmetric part of 
the displacement gradient) in the physical domain does not transform into the sym- 
metric part of the displacement gradient in the reference domair@. Consequently, 
the analysis for affine meshes cannot be simply reproduced for curvilinear ones, 
and new interpolation operators have to be carefully drafted. We have managed to 
prove only an asymptotic stability for the curvilinear meshes. 

^The same problem is encountered in the case of complex stretchings introduced by Perfectly 
Matched Layers, see e.g. [13] . 
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(a) Adaptive /i-refincments, p = 0. (b) Adaptive ft-refinements, p — 1. 




(c) Adaptive /i-refinements, p = 2. (d) Adaptive /tp-refinements. 



Figure 9. Comparison of FE error for adaptively refined curvi- 
linear meshes with the best approximation error. 

Presented numerical results on /i-adaptivity go beyond our analysis as we use 
1-irregular meshes with hanging nodes supported by an existing hp software. 

Finally, the paper presents only a two-dimensional result. We continue working 
on the 3D case using different ideas and hope to present new results soon. 

Appendix A. Mesh generation 

We assume that the domain f2 is a (curvilinear) polygon, and that it can be 
meshed with a regular family (Th)h of C°-compatible meshes (i.e. f2 = UtgT,, ^' 
for all h) that satisfy the regularity assumptions discussed in the previous section. 

We will outline now shortly how one can generate such meshes in practice. 
Suppose the domain Q has been meshed with a C°-compatible initial mesh 7I n t, 
Q = (J*li Ti, where {X^}™! are curved triangles of class C 1,1 . We denote by 
{Gi,-- - ,G m } the mappings from f to {T u ■ ■ ■ ,T m }. Then G { £ C 14 (f) and 
G^ 1 G C 1 (Ti) for any 1 < i < m. For examples of techniques to generate an initial 
mesh satisfying the assumptions above, see |12j . 

Lemma A.l. Let T = G T (T) be a closed triangle in R 2 with G T S C l ^{f). Let 
& > be a positive constant. For any h > 0, we denote by T- h any triangle contained 
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Figure 10. Convergence history for adaptive h- and hp- 
refinements on curvilinear meshes. 



in T such that the diameter of 'Ti is h, and hj p < a where p is the diameter of the 
sphere inscribed in . Let T 3 x — > H x = Bic + b be an affine mapping from T 
onto f- h . Let p be the centroid off. We put B = D(G T oH)(p), b = (G T oH)(p), 
and *(x) = (G T o H)(x) - B(x -p)-b. 
Then, we have 

(sup\\D*(x)\\)\\B-i\\<Ch, 
where C is a constant independent of h, and Ti. 

Proof. Obviously, B = DG T (H(p))B. So B^ 1 = B' 1 (L>G T 0ff (p))) -1 • We have 

D*(x) = (DG T (H(x)) - DG T {H(p)))B. 

Since Gt is a C 1 -diffcomorphism from T onto T, we have sup Xg y, || (Z3Gt(x)) _1 || < 
oo. With G T € C M (T), we also have 

pG T (if(x)) - DG T (H(p))\\ < M\\H(±) - H(p)\\ for any x G f 

where M is the Lipschitz constant for all first order derivatives of Gt on T. 

With hip < a, we obtain that IIBII ■ < -, for some c > 0. With h 

a 

denoting the diameter of T~ h , we have ||-ff(x) — -ff(p)|| < h for any x <E TV. The 
definition of $ implies then 

(sup llD^lDpr 1 !! < |M(8U P UpGrCy))- 1 !!)^. 

xGT 17 yef 

Setting G = -^M(sup- f || (T^G-r (y) ) _ 1 1 1 ) finishes the proof. □ 
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Let Ti be now multiple copies of the reference triangle corresponding to the initial 
mesh, Ti = Gi{Ti), i = 1, . . . , m. 

Definition A. 2. (Regular triangular meshes in the reference space) A family of 
triangulations (7~ of reference triangles Tj is said to be regular provided two 
conditions are satisfied: 

(i) Partitions of edges of Tj mapped into the same edge in the physical space 
are identical. 

(ii) sup^ supj, g -j- h/p < oo, where h and p are the outer and inner diameters 
of f. 

Obviously, uniform refinements of reference triangle^ are regular. A number 
of adaptive refinement algorithms produces regular meshes as well. To this class 
belong e.g. Rivara's algorithm (bisection by the longest edge), Arnold's algorithm 
(bisection by the newest edge), the Delaunay triangulation (see [8]). 

Using Lemma [A. II and the fact that Gi is C 1 -diffeomorphism from T to Tj, 1 < 
i < m, we easily conclude that any regular refinements in the reference space 
produce curvilinear meshes that satisfy our mesh regularity assumptions. 

Appendix B. Properties of Sobolev spaces on curved and reference 

TRIANGLES 

Lemma B.l. Let T be a curved triangle. For any lo £ if^TjR 2 ), we define w(x) 
on T by 

DGr(x) . . - 
" (x) = det(ZX? T ( X ))" (x) XG 
Then w(x) G H(T). Divergence transforms by the classical Piola's rule: 

diwj(x) — (det(DG , T(x)))~ 1 (ii'yxt2;(x) 
for x G T almost everywhere. 

Proof. Notice that w(x) = det(L>GT(x))(DG T (x))" 1 cj(x) for any x G f . It is 
straightforward to see that det(-DGT(x))(L>67T(x)) _1 is a matrix whose entries 
contain first order partial derivatives of Gt{x). Notice that 

w(x) = det(L>G T (x))(L>G T (x))- 1 a;(x) 

is the standard pull back mapping from H (div, T) to H (divx, T). So we immediately 

have divw(x) = - — . — — — ; — rrdiv x c2;(x) for x G T almost everywhere. Since Gt is 
v ; det(DG T (x)) w 

a G 1 -diffeomorphism from T to T, we can conclude that w(x) G 7J(div,T). 

Let e G Ai(T). We denote by ((s) : [0, 1] — !> e, the local affine parameterization 

of e. We have then 

2 _ / r,",r/-r„\\\T/ 



J[0,1] 



(dct( J DGr(C(s)))) 2 W (G T (C(s))) 1 [DG T (t(s))- 1 DGrfCW) - 

0,1] 

W (G r (C( S )))||(£>G T (C( S )))- 1 [ J DG T (C( S ))C( S )]||^. 



4 With the same number of divisions for each reference triangle. 



MIXED FEM FOR 2D LINEAR ELASTICITY 



35 



Since Gt is a G^diffeomorphism from T to T, we can conclude that ui\ g f G 
L 2 {df;R 2 ). So £) € H(f). □ 

Lemma B.2. There exist 8 > and C > sucft. £fta£, /or any h < S and T G 7ft,, 
IMU*(T) < C||w|| i2(f) , V W £L 2 (T;M 2 ), 

where - ^j^q "p"J) = w ( x ) ^ or arl ^ x G P. If T is a triangle for any T G 7ft and 

/i > 0, iften £fte above inequality holds for any ft > 0. 

Proof. 

IMUs(T) = jf det( ^ rW)2 [^GT(x)^(x)] T /JG T (x)^(x)dx 
= [det((DG T )~ 1 )u T [DG T ~DGT]tidx 

= / det( J DG T 1 J B T )det(B- 1 )tZ; T J Bj[ J B 7 ; T J DGf J DGTS T 1 ] J B T c2;o!x. 

Since Cft — >■ as ft- — > 0, lim^o sup TeTji sup Ae y \\B^ T DG^DGtB^ 1 — I\\ = 0. 
By LemmaEU we have lim/j_»o sup Te7 - h sup ig j, | det^G^Br) — 1| = 0. 

Since (7~h)h is regular, there is a constant it > such that o"i/<7 2 < cr for any 
P G 7ft, where o"i and o"2 denote the biggest and smallest singular value of the 
corresponding matrix Bt- Then ||Pj|| = ||-Bt|| = o"i and det(i?T) = 0"i • 0"2- So 
||Pj|| • \\B T \\ det( J B 7 ; 1 ) = 0-1/0-2 < a for any ft > and any T G 7ft. 

We can conclude thus that there exist 6 > and G > such that, for any h < S 
and T G 7ft, 

M L 2 {T} <C\\tj\\ L2{t) VwGP 2 (P;R 2 ) 

If, for all ft, P G 7ft are (regular) triangles, the asymptotic argument is not necessary, 
and the above inequality holds for any ft > 0. □ 

Lemma B.3. There exist 6 > and G > such that, for any ft < 5 and T G 7ft, 

where T (~) -tttt = wfx) /or an?/ X G P. If T is a triangle for any T G 7ft and 

det(_DGT(x)) 

ft > 0, iften ifte above inequality holds for any ft > 0. 
Proo/. 

IH| 2 2(f) = ^(det(CG T (x))) 2 W (x) T [i9GT(x)- T ^GT(x)- 1 ] W (x)dx 
= / det(DGT(x))w(x) T [P>G T (x)- T P'GT(x)- 1 ]a;(x)dx 

JT 

= ( det(DG T (x)P T 1 )det(P T )a;(x) T B T T 

JT 

[B^DGT{xy T DG T {xy 1 B T }B T 1 uj(x)dx. 

According to Lemma 17^1 linift^o su PTeT h su Pxef W^t -DGt(x)~ t P>Gt(x) _1 Bt — 
J|| =0. By Lemma l7.9l we havelhm^n sup Tc7 ^ sup f . c j. | det (DGt (x) 5 j 1 ) — 1| = 0. 
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Since (Th)h is regular, there exists a constant a > such that <7\/<J2 < cr 
for any T £ 7h, where cri and 02 denote the biggest and smallest singular value 
of the matrix 5 T . Then || J B T T || = \B T l \\ = a^ 1 and det(S T ) = <Ji • cr 2 . So 
||£?y T || • H-B^H det(Bx) — (Ji/<72 < cr. Consequently, there exist <5i > and G\ > 
such that, for any h < 8\ and T £ 7/j, 

det(Z?G T (x)B- 1 )det(S T )w(x) T B T T [BjZ)GT(x)" T Z?GT(x) _1 J B T ] J B T 1 a;(x) 
< G 2 cj(x) t u;(x) 

for all u £ M 2 ), x £ T. We can conclude that, for any h < 5\ and T £ 7ft, 

||w|| i2(f) <Ci|MU2(t) V W £H 1 (T;R 2 ) 
According to Lemma |B.1| we have 

||divxw|| 2 2(f) = y (det(L>G T (x))) 2 (divw(x)) 2 rfx = ^ dct(L>G T (x))(divw(x)) 2 dx 

= / dct(L»G T (x) J B T 1 )dct( J B T )(diva;(x)) 2 dx. 

JT 

By Lemma 1731 lim/^o sup TeTh sup ig j, | det(DGT(x)i?^ 1 ) — 1| = 0. Obviously, 
dct( Bt) < h^- There must exist then <5 2 > and G 2 > such that, for any h < <5 2 
and T £ 7fc, 

det(Z?G T (x) J B 7 ; 1 )det( J BT)(diva;(x)) 2 < G|/i|(diva;(x)) 2 £ iJ^TjM 2 )^ £ T. 
We conclude that, for any h < 82 and T gTh, 

\\div A U)\\ L2if) < G 2 /i T ||divw|| L 2 (T) £ i/ 1 (T; R 2 ). 

We take now an arbitrary e £ Ai(T). We denote by £(s) : [0, 1] — > e the local afhnc 
parameterization of e. We have then 

M\l H e) = / (^(C( S ))) T ^(C( 5 ))IIC( S )II^ 

-'[04] 

= / (dctiDGTiasMMGHC^ViDGTias^DGTias))- 1 } 

u(G T (as)))U(s)\\ds 
= I (det(DG T (C(s)))B^) 2 (det(B T )) 2 Lj(G T (C(s))) T B^ T 

J [0,1] 

[^^G T (C(s))- T ^G T (C(s))- 1 ST]i? T ^(G T (C( S )))||C( S )||rf.s. 
According to Lemma 17. 81 we have that 

lim sup sup sup \\B^DG T {C{s))~ T DGriCis^Br - I\\ = 0. 

h ^°T£T h eSTse [0,1] 

By LemmaEl lim^o sup T6Tft sup eeT sup se[0J] |(det(L»G T (C(s))) J B T 1 ) 2 - 1| = 0. 
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Consider again the singular values of Bt, ci > <?2- Then ||i3 T T || = 1 1 S T 1 1 1 = 
cr^ 1 , and det(B T ) = en ■ ff 2 . So ||B T T || • || J B T 1 ||(det(B T )) 2 = (cti) 2 < h 2 T . Conse- 
quently, there exist S3 > and G3 > such that, for any h < S3 and T £ 7ft, 

(dctpG T (C(s)))i? T 1 ) 2 (det( J B T )) 2 W (GT(C(s))) T BT T 
[Sf J DGT(C(s))" T ^G T (C(s))" 1 S T ]B T 1 ^(GT(C(s))) 
<G 3 /iTw(x) T a;(x) Vw £ ff^TjR 2 )^ e Ai(T),sG [0,1] 

and 

||C(*)[| =||^ 1 [5 T £»G T (C( S ))- 1 ] • [DGr(C(5))C(«)]|| 

<G 3 ^II^MC(s))C(s)ll Ve £ A x (T),a £ [0,1]. 
We can conclude that, for any h < S3 and T £ 7ft, 

W\h W < Cfh T I c(G T (C( S ))) T c(G T (C( S )))pG T (C( S ))C( S )||d S 

J [0,1] 

= ^rlMli» W - 

Obviously, sup^ sup Ter(i hr < 00. Since lj £ i/ 1 (T;R 2 ), we can use the Trace 
Theorem to conclude that there exist S > and G > such that, for any h < S 
and T € 7ft, 

It is easy to see that, if T is a (regular) triangle for any T £ 7ft and h > 0, then 
the inequality above holds for any h > 0. □ 

Lemma B.4. There exist 5 > and G > smc/i i/iai, /or any h < S and T £ 7ft, 
\\curhj\\ 2 L2{T) < Ch T 2 \\curkoj\\ 2 L2(f) Vcj £ iJ^TjM 2 ) 

where T " ^ " r = cj(x) /or any x £ T. If T is a triangle for any T £ 7ft and h > 0, 
det (-Bt) 

iften £/ie above inequality holds for any h > 0. 
Proof. We have, 

CUrlCj(x) = d^ (cUrL2>)(A) = det(i? T ^G T (x)) i3TCUrl ^ ( " )(JGT(jc))T 
- (det (_By 1 )) 2 .^curi^^sTp^^^-^T 



det(B^ i J DG T (x)) 



and 



|curl^|| 2 2(T) = / J^i n^^a PGr^^^rH.^x)) 1 
v ' Jt (det(_B T Z?Gt(x))) 2 

[Bf B T ] curl x w (X) (L>G T (x) B T 1 ) T dx 

(^^t 1 )) 3 (Z9G T (x) J B 7 : 1 )i?T(curl i ^(x)) 

f (det^^GrW)) 2 

[B? B T ] curlxd; (x) B^(DG T {^)B^ 1 ) T d±. 

Since Cft — )■ as ft. — > 0, lim^o sup Te7 ^ sup ieT ||£)Gr(x)Sy 1 — J|| = 0. 
By Lemma I7T91 linift-x) su PTer h su Pief I (det (-B^ 1 Z?Gt(x))) 2 — 1 =0. 



T 
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Since (Tk)h is regular, there exists a constant a > such that o~x/o~2 < a for 
any T € 7ft, with u\ > o"2 denoting the singular values of matrix Bt- Then 
||B/|| = H^xll = ax and det(B^) = <rf 1 • a^ 1 . So || 2 • H-BTfCdet^ 1 )) 3 = 

ci/cf < o"/ "! — c '^t 2 f° r some constant c' > 0. 

Consequently, there exist Si > and C > such that, for any h < 8\ and T G 7ft, 

, ( det (V)) 3 (^G T (x)B 7 ; 1 )B T (curl^(x)) T 

f (detOB^ZXMx))) 2 V V ' T ' V V 

[Sj S T ] curUw (x) (L»G T (x) B T 1 ) T dx 
<Ch^ 2 / ||(curl x w(x)) T (curl x w(x))||dx. 

Again, it is easy to see that, if T is a (regular) triangle for any T e % and /i > 0, 
then the above inequality holds for any h > 0. This finishes the proof. □ 

Lemma B.5. There exist S > and C > smc/i £/ia£, /or any h < 5 and T G Th, 

IM&ccu*,*) <c(IMIi 2( T) + ^Tll«^MIi 2( T)) vweff^K 2 ) 

where - T " r = cj(x) /or an?/ x G T. If T is a triangle for any T G 7ft and /i > 0, 
det (Bj>) 

then the above inequality holds for any h > 0. 
Proof. We have, 



and 

|2 



MI^T) 



: / (det(B T )) 2 w(x) T J B T T J B T 1 w(x)dx 

= / det(B T )det( J B Tj DG T (x)" 1 )w(x) T J B T T B T 1 w(x)dx. 
Jt 

By Lemma l7Jl limft_>o sup Te7 - h su Pxgt I det^T-DG^x) -1 ) — 1| =0. 

Since (7~h)h is regular, there exists a constant a > such that o~i/o"2 < 0" for 
any T G 7ft, where o~i > o"2 are the singular values of matrix i?T- Then ||_B/ T || = 
\\B^\\ = o-^ 1 and det(B T ) = <*x • o- 2 . So, ||S 7 ; T || • HB^ 1 !! det(B T ) = 01/0-2 < o\ 
Consequently, there exist 61 > and Gi > such that, for any h < 8\ and T G 7ft, 

dct( J BT)det(B T Z?GT(x)" 1 V(x) T J B T T B T 1 a;(x) < G 2 w(x) t cj(x) 

for any cj G R 2 ), x G T. We can conclude that, for any /i < 0/ and T G 7ft, 

IHI i2(T) < Cx\\uj\\ lHt) ,Vuj G ff^R 2 ). 

At the same time, 

curl x a)(x) = det(_BT)S/ 1 curl x w(x) 

= det(B T )B T 1 curl x w(x)(£)GT(x))" T det(LiGT(x)) 

= det(B T )B T 1 curl x w(x)B T T (BfL>G T (x)" T )dct(Z)GT(x)) 

= det (B T ) B T 1 curl x lj (x) B t t (DG t (x) " 1 B T ) T det (DG T (x) ) . 
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and 

||curl*w|£ a(f) =^det(B T ) 2 detpGT(i)) 2 ( J DG T (i)" 1 S T ) J B T 1 (curl xW (x)) T 
B T T B T 1 curl x a; (x) B T T (DG T (x) " 1 #t ) T di 

= [ det(B T ) 2 det(DGT(x))(L»GT(x)" 1 J B T )B 7 ; 1 (curl x w(x)) T 
B T T B T 1 cur l x w (x) B T T ( DG T (x) " 1 B T ) T dx 

= / det(B T ) 3 det(B T 1 £)G T (x))(Z)G T (x)- 1 J B T ) J B T 1 (curl x a;(x)) T 

S T T B T 1 cur l x w (x) B T T ( DG T (x) " 1 B T ) T dx . 
According to Lemma lim^_^ sup TgTh sup ig j, || DGt(x) _1 .Bt — I\\ = 0. 
By Lemma \T§[ lim^o su PTer h su PxeT I det (Sy 1 Z3GT(x)) — 1| = 0. 
Since (7h)/i is regular, det(£?T) 3 ||-BT T || 2 ' II-^t 1 !! 2 — for some constant a > 0. 
There exist thus 62 > and G2 > such that, for any h < 5 2 and T € Th, 

det (B T ) 3 det (B T 1 £>G T (x) ) {DG T (x) " 1 S T ) 5 T 1 (curl x cj (x) ) T 

B T T J B T 1 curl x w(x)B T T (L>G T (x)" 1 B T ) T dx 

<G|/(,|||(curl x w(x)) T curl x w(x)|| Vw £ _ff 1 (T; R 2 ), x 6 T. 

We conclude that, for any h < 62 and T £ 7h, 

||curl x w|| i2(f) < G 2 h T 1 1 curLv 1 1 L 2 (T) ,Vw e ff^TjR 2 ). 

Again, it is easy to see that, if T is a triangle for any T £ Th and ft. > 0, then the 
inequality above holds for any h > 0. This ends the proof. □ 



u = u o G T 1 or u — u o Gt 



Lemma B.6. Let T be a curved triangle. Let w(x) fee defined on T and u(x) fee 
defined on T and 



where Gt is the element map. Then u(x) £ L 2 (T) if and only if u(x) £ L 2 {T). 

Proof. This is an immediate consequence of the fact that Gt is a G^diffeomorphism. 

□ 

Appendix C. Proofs of results from Section 8 
Proof of Lemma [ 



Proof. We put r = f(T). We assume {£i(x),--- ,£(x) ir } is a basis for VfA 2 (f). 
Then 



n? T u(x(x)) 



1 lr 
det(£>G T (x))§ a ^ i(i) 



and 



^(x(x))= det( ^ T( , )) gM(x) 

Coefficients (ai, • • • , a^) T and ■ • • , /3z,.) T are obtained by solving the following 
two linear systems, 

Ax (ai,--- ,aO T =bi and A 2 (/3i, • ■ • , A r ) T = b 2 
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with 

&(x(x))fr(x(x))rfx 
l T det(DG r (x(x))) J T 



Ct(x(x))| 3 -(x(x))tix /" £j(x(x))dx 



(j4a) « - i T (det(DGr(x(x))))=» ^ ~ 7 X det(£>G T (x(x))) " (x) 
for 1 < i, j < l r . By pulling back to T we obtain, for any 1 < i, j < l r , 



f &(x)£ i (x)dx (bi^ = det(fl T ) / det(S T 1 ^G T (x)) U (x)4(x(x))dx 

Since det(£>T) is a non-zero constant, (det (Bt) A 2 ■ • ■ , Pi r ) T = det(i?T)b2- 
So we can redefine A 2 and b2 in the following way. 

{M)a = J d ^GTm ^ j{i)dk (b2)j = dct{BT) L u M k »&w dk 

According to Lemma \THl limh-x) sup rg7 - fe ||Ai — A 2 \\ = 0. And, for any e > 0, 

lim sup ||bi -b 2 || 2 <(edet(B T )) 2 / u 2 (x(x))dx 
h^a TeTh Jf 

=e 2 dct(B T ) [ dct( J B T L>G T (x(x))- 1 )w 2 (x)dx 
Jt 

<4e 2 det(Br)||«||| 3(T) . 
The last inequality holds when h is small enough. This implies that 



lim sup ||bi -b 2 || < 2e y / dct(B T )\\u\\ L 2, T) . 
h-M TeTh 

So, for any e > 0, there exists S > such that, for any h < 5 and T £ Tk, we have 
||(ai-/3i,-- - ,a,„ -AJII < 3s^det(B T )\\u\\ L 2 iT) . 

We have then 

l|n 2 T u - p, t u||| 2(t) = jf (det( ^ (x))))2 £ £(*M)(a, - A) 2 rfx 

= /, detpG T (x)) | l?( " )(Qi - A) ^ 

<ce 2 |M|| 2(T) / det^X^x)- 1 )]^*)** 
Jf i=i 
Here c is a positive constant which depends on Z r only. By Lemma 17.91 there exists 
M > such that, for any h small enough and any T € 7ft, 

[ det(B T £>G T (x)" 1 ) Ve 2 (x)dx < M 2 
J T i=l 
We can conclude that, for any e > 0, there exists <5 > such that, for any h < S 
and T e 7/ l7 

||IIp iT « - Pf,Tu|U 2(T) < E||«|| £a(T) ,Vu e i 2 (T) 
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Here P?,t is the standard I/ 2 -projection onto V?A 2 (T). 
Proof of Lemma EU 



□ 



Proof. Obviously, d>(x) = det(DGT(x))(DG T (x.))~ 1 u(G T (x)). Using Lemma EZQ 
we can conclude that w(x) G H(T). 

By pulling back to T and using the definition of Vf+ih}(T), we can see that (I9.4[) 
is the same as (|9.7[) . and (|9.3j) is the same as (|9.6[) . Thus, we only need to show 
that (|9.5j) is the same as (|9.8j) . Since 7h is C°-compatible, then Gt(C( s )) = x e(s) 
for any s G [0, 1] where £ : [0, 1] — > e is an affine local parameterization of e. 
We have then 



[0,1] 



[w(xe(s)) • n(x e (s))]r}(s)||x e (s)||ds 



[0,1] 



[0,1] 



[0,1] 



[w(G T (C(a))) • n(Gr(C(«)))]^(*)|pGr(C(«))C(»)ll*» 



£>G T (CW)w(C(a)) ( J DG T (C( S )))- T n(C( S )) 



det(£>G T (C( S ))) || (DGt(C00))- t A(CW) II 
<&(£(*)) A(C(*)) 



^(a)||UGr(C(a))C(*)[[ds 



det(UGr(C(*))) ||(^G T (C(s)))- T n(C( S 



r?( S )||£>G T (C(s))C(s)||rf S 



Notice that £(s) = ci where t is a unit tangent vector along e, and c is a nonzero 

constant, n(C(s)) - (t 2 ,-ii) T , and (DG T (C(s)))- T = ., n< ^,„ with 

det(DG T (C(s))) 



(DG T ) 22 (C( S )) -(DG T )ai(C(«)) 
-(£>G T ) 12 (C(s)) (DGt)h(CW) 



Therefore, we have 

«(C(a)) 



[0,1] 



[0,1] 



det(£>G T (C( S ))) ||pG T (C(.s)))- T n(C( S ))| 
[<&(£(«))• fi(C(*))]»)(»)dfl 



f7(*)||DG T (C(*))C(*)l|da 



We conclude that (|9.5|) is equivalent with (|9.8p . This finishes the proof. 

Proof of Lemma 19.101 
Proof. By Lemma iB.li we have w(x) G H{T), and 

divw(x(x)) = - — - — — — — — — divx(l)(x) 



□ 



det(DG T (x)) 

for x £ T almost everywhere, provided we define w(x) on T by 

DGt(x) 



;(x(x)) 



det (£)G T (x)) 



w(x) 



for any x G T. 
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divll~ +1 T w(x(x)) = div 



Using Definition 19. 1[ Definition 19.21 Lemma HQ] Lemma T9.91 and Lemma 10 in 
[19], it is easy to see that 

H? T div W (x(x)) = dct(g g r(A)) ^,fdiv ft a(x) = ^-^i-^div^ ^(x) 

L>G T (x) , 1 , , 

. det(DG T (x)) n ^^" (x) J = det(DG r (x)) dlV ^^^ (x) 

We have thus II? T divw = divlll , x T w. □ 

Proof of Lemma 19.111 

Proof. According to Lemma IB.2I and Lemma IB.31 there exist 5 > and C\ > 
such that, for any h < 8 and T G 7~k, 

||IIf,Tw|U»{T) < Ci ||n^|| i2(t) G i 2 (T; K 2 ) 
ll^ll 2 ff( div„T) + IMI^af) ^ Cill^ll^m e ^(TjK 2 ) 

By Lemma MM ^ £ iT^R 2 ) for an Y w e iT^R 2 ). 

The definition of operator II- f implies that there exists a constant C 2 > such 
that 

j{ny^(i)) T ny^)dk < c 2 (\\uf H{d . v ^ f) + ||^|| 2 2(e(f)) ),v^(x) g [H\f )] 2 . 

It is easy to see that for affine meshes the above inequality holds for any h > 0. 
This finishes the proof. □ 



Proof of Theorem 19.151 

Proof. For any h > and any T G Th, we define a linear isomorphism At from 

H 1 (f ;M 2 ) to i^fTjR 2 ) by (A T w)(x(x)) = f ■ It is easy to sec that A T is a 

det^-Dr) 

linear isomorphism from Vf. +2 A°(T; R 2 ) to P f+2 A°(T; R 2 ). 

We define an operator £ T : R 2 ) — > P f+2 A°{f ; R 2 ) by W T (A T w) = 

j4r(£rw). Obviously, Wr is well-defined if and only if Et is well-defined. We 
denote by {£1, • ■ • , £;.} a particular basis of Pf + 2A°(T; R 2 ). 

According to Lemma fB. 41 and Lemma fB.SI it is sufficient to show that there exist 
S > and C\ > such that, for any h < 8 and T & Th, Et is well-defined, and 



(zi, ■ ■ ■ ,zi?)\\ < Ci\\lo\\ h1 ^ for any Cu. Here J2k=i z k(,k = Et 



to. 



According to the definition of Wt, Et can be defined by relations 
(C.l) f dW(A T E T 0J - A T u){x)$(x(x))dx = V</> G Vf (T )(T)/M. 



(C.2) 

J {(A T E T u>)(x) - (A T c2;)(x)) T ^GT(x(x))- T h l (x(x),i j;(f) )dx = l<i<k? 
(C.3) 

[(A TJ B T w-A T w)(x e (s))-n(x e (s))]r7(s)||x e (s)||ds = V?) G P f(e) ([0, l]),Ve 6 A X (T) 



[04] 
(C.4) 

/ [(yl TJ B T w-A T w)(x e (s))-t(x e (s))]^(s)||x e (s)||ds = V77 G 7W[0, l]),Ve G A X (T) 

J [0,1] 
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(C.5) E T u = at all vertices of T 

Denote: 



hi bi2 

&21 &22 



B T ,J = det(J3 T ), 



Ul 
U-2 



Wi 
W-2 



By pulling back to T, Et can be defined by relations 
(C.6) 

/ J' 1 [(b 11 {DG T )22 - 6aipGr)ia)fii,i + {b l2 {DG T )22 - b 2 2(DG T ) 12 )u2,i 
Jf 

+ (b 2 i(DG T )u - bu(DG T )2i)ui,2 + (b 2 2(DG T )n - b 12 (DG T )2i)u2,2}ip{x)dii 

= / J" 1 [(6i 1 (^G T )22 - M-DG T )i 2 )u>i,i + (bi2(DG T )22 - M£G T )i 2 )u>2,i 
Jf 

+ (b2i(DG T )u - &ii(DG T )2i)t&i,2 + (b22(DG T )u ~ b 12 (DG T )2i)u> 2 ,2]4>(^)dx, 
V^GP f(T) (f)/f 

(C.7) /(^ T w(x)) T Bji?G'T(x) _T h l (x,i)det(B T 1 i?GT(x))dx 

= /(w(x)) T SfZ?GT(x) _T h l (x^)det( J B T 1 J DG T (x))'ix 1 < i < fe? 



(C.8) / [5 t £tw(C(s)) • (DG r (C(a))" A(C00))tf00 det(S- 1 £»G T (C(s))^ 



[B T w(C(a)) • (DG T (C(s))- 1 n(C(s)))]r}( S )det(^ i DG T (C( S ))d S 

'[0,1] 

Vi)€Pf W ([0,l]),Ve€Ai(T) 



(C.9) 



[Bt^tw(C(s)) • (^G T (C(s))C(s))]'7(s)det(^ 1 )d S 

[0,1] 

[BrtS(C(*)) • (DG T (C(s))a S m( s )det(B^)d s Vr) € 7> f(e) ([0, 1]), Ve € A X (T) 

[0,1] 

(CIO) £ T w = at all vertices of f 

It is easy to see that (IC6|) comes from (IC.1[) . (|C.7I) comes from (|C.2I) . (|C.9[) comes 
from ([(13)1 . and (jC.10|) comes from ([HI)]) . And (|U15|) can be got from by 
using the fact that ||x e (s)|| = \\DG T (C(s))i\\ = c||DG T (C(s))" T n|| for some non- 
zero constant c, which comes from direct calculation. 

Notice that vector ((s) is constant tangent vector along each edge of T. Set 

C(s) 1 B f B T ({s) 

Obviously, 6/0. 
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Perform now the operation: b x [ (|C.9j) -ax (1C.8|) ]. We have, 



(C.ll) / [3ru(£(a)) ■ F T (a)]fj(a)da = / [u>(C(a)) ■ F T (a)]r}(a)da 

J[0,1] J[o,i] 

V77GP f(e) ([0,l]),Vee A^T) 

where 

F T (a) = det(B^)[B^(DG T (as)))C(a) 

defc(Br)||C(*)|| 



det(Br 1 Z>Gr(C(«)))(n 1 B^B T Q{a))B^ (DG T (C(s)))~ ' n)] 



asVB^B T C(s) 

Then the definition of operator Et can be rewritten by using conditions (|C.6|) . (jC.7|) . 
(1081) . (ICHl) . and (IClUt . 

Using the fact that n±£(s), Lemmas l7.7l7.8|7.9l and the assumption that (Th)h 
is regular, we obtain 

C(s) 

lim sup sup sup \\Ft{s)-—. 1|| = lim sup sup sup ||F r (s) -n|| = 

h^0 TeTheeAl{T)se[0A ] ||C(s)ll h ^°TeT h e€A 1 (T)ae[0,l] 



Consequently, 



lim sup sup sup ||Ft(s) — t|| = 

h^0 TeTh eeAi(T) sG[0,l] 



We denote now by E(T,r) the matrix corresponding to the left-hand side of con- 
ditions (fCT6 |) .(fCT7 1) .(fCT8 |l . ([CTTT) . ([CTT0|) . a particular basis of V,, +2 A°(f ; R 2 ) (the 
solution space), some basis of V? (t) (T)/^, and some basis of V?( e )({0, 1]) for each 
e £ Ai(T). We denote by {&,••• ,6,,} a basis for P f+2 A Q {f-R 2 ). Finally, we 
denote by C(f) the matrix corresponding to the left-hand side of conditions (|5 . 14[) . 
()5.15|) . ()5.16p . and (|5.17|) . and the same bases as above. 
Using the fact that 

lim sup sup sup ||-FV(s) — 1|| = 

h~>0 TeTh eeAi(T) s£[0,l] 

and Lemmas I7.7|7.8[ l7.9[ we conclude that, for any t £ [0, 1], 

lim sup \\E(T,r) - C(r)\\ = 
h^o TeTh 

Then, for any given t £ [0, 1], and any given f with non-singular C(f), the matrix 
E(T, f) is non-singular for any T £ Th when h > small enough. Notice that the 
right-hand sides of conditions (|C.6[) . (|C.7|) . (|C.8|) . (|C.lip . and (IC.lOj) are continuous 
linear functionals of w £ i? 1 (T;IR 2 ). We can conclude thus that the operator Wr,t 
is well-defined for any T £ Th with small enough h. 

Since for any t £ [0,1], lim^o su PreTh ll-^C^^) ~ CC^ 1 ) 1 1 = 0j an d the matrix 
C(f) depends only on f, we can conclude that there exists C\ > such that when 
h > small enough, then ||(zi,-- - < C\ \\^\\ H iif\ for any T £ Th- Here 

Finally, it is easy to see that the operator Wt will be well-defined and the 
inequality in the statement of this theorem holds for any h > for affine meshes. 
This finishes the proof. □ 
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